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ABSTRACT 


The Saturnian radiation belts have not received as much attention as the Jovian radiation 
belts because they are not nearly as intense — the famous Saturnian particle rings tend to deplete 
the belts near where their peak would occur. As a result, there has not been a systematic 
development of engineering models of the Saturnian radiation environment for mission design. 
A primary exception is that of Divine (1990). That study used published data from several 
charged particle experiments aboard the Pioneer 11, Voyager 1, and Voyager 2 spacecraft during 
their flybys at Saturn to generate numerical models for the electron and proton radiation belts 
between 2.3 and 13 Saturn radii. The Divine Saturn radiation model described the electron 
distributions at energies between 0.04 and 10 MeV and the proton distributions at energies 
between 0.14 and 80 MeV. The model was intended to predict particle intensity, flux, and 
fluence for the Cassini orbiter. Divine carried out hand calculations using the model but never 
formally developed a computer program that could be used for general mission analyses. This 
report seeks to fill that void by formally developing a FORTRAN version of the model that can 
be used as a computer design tool for missions to Saturn that require estimates of the radiation 
environment around the planet. The results of that effort and the program listings are presented 
here along with comparisons with the original estimates carried out by Divine. In addition. 
Pioneer and Voyager data were scanned in from the original references and compared with the 
FORTRAN model’s predictions. The results were statistically analyzed in a manner consistent 
with Divine’s approach to provide estimates of the ability of the model to reproduce the original 
data. Results of a formal review of the model by a panel of experts are also presented. Their 
recommendations for further tests, analyses, and extensions to the model are discussed 

KEY WORDS: Saturn, Radiation Models, Radiation Belts, Cassini Spacecraft, High Energy 
Electrons, Trapped Particles, Space Radiation 
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EXECUTIVE SUMMAR Y 


The primary objective of the Prometheus Saturnian Radiation Modeling task was to 
develop engineering models of the radiation environment near Saturn based on data from Pioneer 
1 1 and from Voyagers 1 and 2 so that Prometheus mission planners could estimate the degree of 
radiation exposure for missions to the Saturnian system. This was accomplished by collecting 
and reviewing available high-energy electron and proton data relevant to the Saturnian system as 
originally done by Divine [1], In that original study. Divine constructed models of the electron 
and proton fluxes as functions of pitch angle (a), magnetic field at the position (B), and magnetic 
field line passing through the location (L-shell). The current effort entailed encoding Divine’s 
original algorithms and developing methods for estimating the Saturnian magnetic field 
components (B and L) needed to estimate the fluxes along a spacecraft trajectory. While 
originally used as an engineering tool to model Cassini spacecraft doses during its design, the 
current model implementation (called SATRAD) was used to estimate the doses for Prometheus 
missions to Saturn and will be used to compare with the Cassini measurements. To validate the 
SATRAD model, its fluxes have been compared by the authors to published Pioneer and 
Voyager data in the same manner as Divine. Those comparisons have been used to estimate the 
“accuracy” of the SATRAD model by calculating an “error range.” Within the limits of the 
current study, the original estimates of Divine have been validated (e g., to a factor of 2-3) and 
are well within current expectations for radiation modeling tools. 

To independently validate the model and at the request of the Prometheus Program, a 
review was carried out by experts in the area of Saturnian magnetospheric modeling on 18 
March 2005 to determine how well the model performs. While the Review Board found the 
model to be a useful baseline, the members provided a range of suggestions on how the model 
could be improved in the future. The Review Board comments on the model can be roughly 
divided into 4 groups: 

1) Pioneer and Voyager particle data: 

a) Need to evaluate using data in the Planetary Data System (PDS). 

b) Need to re-evaluate detector background sources and calibrations in the light of 
current understanding of the instruments. 

c) Pioneer and Voyager data still of value because they cover energy ranges and 
locations not covered by Cassini. 

d) Potentially important local time (relative to Sun) and real-time “weather” 
variations are difficult to separate using only the Pioneer and Voyager data. 

e) Species identification requires more careful consideration — ions versus electrons; 
protons versus oxygen ions; etc. 

f) Some spurious features in model. 

g) Pitch angle distributions need to be reconsidered and a more flexible fitting 
scheme employed. 

2) Cassini data: 

a) Need to incorporate latest Cassini data when it becomes available. Unfortunately, 
only data from the Saturn Orbit Insertion (SOI) orbit may be useful for critical 
pitch angle information as instruments may now be “stuck” in one position. 
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b) Cassini found an inner belt (see Figure 14) in a region that was previously 
unexplored, and thus not included in the model. 

3) Magnetic Field: 

a) Uncertainties in magnetic field make field line tracing difficult as the field 
approaches 0 at ~10 Saturn radii (Rs) — magnetic field models of Saturn are still 
being developed. 

b) There are potentially important local-time (relative to Sun) and real-time (“space 
weather”) variations that are difficult to separate in the data. 

c) Cassini found a 6 minute slowdown in Saturn rotation and a highly variable 
magnetic field in comparison to Pioneer and Voyager observations. 

d) Future magnetic field models will need to include effects of magnetopause and 
ring current in magnetic field model in outer magnetosphere (~10 Rs). 

4) Overall model package: 

a) Need to standardize and have peer reviewed in literature. 

b) Recommend making code/engineering model available for Cassini and other 
programs outside JPL. 

By way of summary, the Review Board found several areas that should be considered in 
both using the model and in updating it. These range from calibration/background issues with 
the original Pioneer and Voyager data to very real concerns over the apparent variability of the 
Saturnian magnetosphere. The reviewers stated that Cassini data should allow improvements in 
the SATRAD magnetic field model and flux estimates and should be actively pursued. Finally, 
the reviewers requested copies of the model for comparisons with the Cassini data and 
recommended that the model be standardized, made more readily available, and formally 
reviewed in the literature as soon as possible. 
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SATURN RADIATION (SATRAD) MODEL 


H. B. Garrett, a J. M. Ratliff, a and R. W. Evans b 

“Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, 
California 91109 

b Gibbel Corp., 2550 Honolulu Blvd., Montrose, California 91020 

INTRODUCTION 

The SATRAD computer model of the Saturnian radiation environment is based 
on earlier work by Divine [1], who examined data from Pioneer 1 1 and from Voyagers 1 
and 2, and used those data to develop an engineering model of the electron and proton 
fluxes over limited energy ranges between ~2 to 10 R s (1 R s = 1 Saturnian radius). In the 
present work, the “Divine model” has been used as a baseline for further development of 
the model and for implementing it into a formal computer code that can be used as a 
numerical tool for engineering calculations of the high-energy particle environment. The 
SATRAD model consists of two components: the original Divine formulation, which 
gives the electron and proton particle fluxes in terms of the so-called B-L coordinate 
system, and a second component which provides the magnetic field model used in 
determining the B-L coordinates necessary for evaluating the model fluxes. To 
demonstrate potential applications of the tool, Fig. 1 illustrates a mission dose/depth 
analysis prepared in support of a Prometheus Saturn mission study. Fig. 2 compares 
radiation estimates for Saturn with similar orbits at Jupiter. To evaluate the predictions 
of the tool, Pioneer and Voyager in situ data at Saturn have been digitized from original 
data plots for comparison with the model. Because much of those data are in count rates, 
count-rate-to-flux conversions were also developed. As an illustration of the process. 
Fig. 3 presents comparisons between Voyager data and the model. The agreement is 
quite good for these energies but other energy channels do not compare as well, implying 
that there are uncertainties in the Saturnian data (see later discussion and appendices). To 
independently evaluate the models and at the request of the Prometheus Project, a Review 
Board was asked to review the different model components, the validity of the data used 
to develop the model, and the ability of the model to predict the actual environment. This 
report details the materials presented to the Review Board and the Review Board findings 
relative to the overall SATRAD model. Finally, as requested by the Review Board, the 
SATRAD FORTRAN model listings are provided for use by the space science 
community. 
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SATRAD Dose 



1.E+00 1 E+01 1.E+02 1 E+03 1.E+04 

Shielding Thickness (mils Al) 

Fig. 1. Estimate of the dose for a potential Prometheus Saturn Moon Mission. The 
full Saturnian magnetic field was used in tracing the magnetic field lines. 


4-Pi Spherical Shell Dose Depth Curves 
1 Week in Selected Circular Orbit 



Aluminum Shell Thickness (mils) 


Fig. 2. Jovian and Saturnian radiation environments for a 1 week exposure. 
Distances are at 2.55 R, 5.95 R, and 9.47 R where R is one planetary radius (Rj = Jovian 
radii, R s = Saturnian radii). Distances correspond to orbits of Amalthea, lo, and Europa 
at Jupiter. Dose in rad(Si) for 4n aluminum spherical shell. 
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Fig. 3. Comparisons between the Saturn model and data from Voyager 2 fly-bys for 
(A) 10<E<20 MeV electron integral flux and (B) 48<E<63 MeV proton integral flux. 
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SA TURN RADI A TION MODEL DE VEL OPMENT 


THE DIVINE RADIATION MODEL 

As a means for describing the development of the SATRAD model, the original 
particle model as discussed in Ref. 1 (N. Divine, JPL IOM-52 17-90-029, Numerical 
Models for Electron and Proton Distributions in Saturn's Radiation Belts, 28 February 
1990) will first be reviewed. Specifically, Ref. 1 is reviewed from the standpoint of 
Divine’s original Saturn radiation modeling process and his methods for estimating 
residuals between the observations and the model predictions. To start. Divine’s Saturn 
radiation model assumes these independent variables: 

B = Magnetic Field strength (gauss) 

L = L-shell [2.3<L<13.0] 

a = pitch angle 

E = Energy (MeV) [0.04<E C <10 MeV; 0.14<E P <80 MeV] 

The following dependent variables are then estimated: 

I = Integral Intensity (cm ' 2 s ' sr 1 ) 
i = Differential Intensity (cm ' 2 s ' 1 sr ' 1 MeV 1 ) 

The algebraic form of the model to predict these variables becomes, in terms of the 
dependence of I on E, 


log I = A 0 -A l log£ + 




l0£ 


E~ + E: 


1 + E 


+ 


2 J 


l0£ 


£ + El 


1 + £ 3 


( 1 ) 


3 J 


and 


dl / , A, -A. A, -A, 

= — A, H H — T 

dE E[ l + (E 2 /E)~ 1 + (EJE) 


Dependence of I and i on a and B is given by: 

A n =y( Si na ) 2 ( a n-b„)+b„ 


( 2 ) 


( 3 ) 


where 

n= 0, 1, 2, and 3 

B c = Magnetic Field strength at Magnetic Equator (gauss) 
a„, b n , E 2 , and E 3 are tabulated at selected values of L — they are to 
be linearly interpolated between the specified L values 
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Values for the various constants are presented in Tables 1 and 2 as taken from Divine 
(Ref. 1). The data and companion references that these parameters were derived from are 
listed in Appendices I, II, and III. Sample output from the model, namely the 
omnidirectional electron and proton fluxes at the equator and at high latitudes, are plotted 
in Figure 4. Also listed is the Saturnian magnetic field model assumed by Divine. 

Table 1. Parameter values for the Saturn electron model (Divine, Ref. 1) 


L a n a , a , a , b „ b , b „ b , E „ E , 

2.30 

4.74 

0.00 

0.00 

6.47 

4.74 

0.00 

0.00 

6.47 

0.60 

1.50 

2.33 

5.02 

0.00 

0.00 

4.01 

5.02 

0.00 

0.00 

4.01 

0.60 

1.50 

2.34 

4.91 

0.00 

0.00 

4.86 

4.91 

0.00 

0.00 

4.86 

0.60 

1.50 

2.36 

4.98 

0.00 

0.00 

4.41 

4.98 

0.00 

0.00 

4.41 

0.60 

1.50 

2.38 

5.00 

0.00 

0.00 

5.37 

5.00 

0.00 

0.00 

5.37 

0.60 

1.50 

2.46 

5.49 

0.00 

0.23 

3.00 

5.49 

0.00 

0.23 

3.00 

0.60 

2.00 

2.56 

5.56 

0.00 

0.34 

3.00 

5.56 

0.00 

0.34 

3.00 

0.60 

2.00 

2.70 

5.51 

0.00 

0.27 

4.48 

5.51 

0.00 

0.27 

4.48 

0.60 

2.00 

2.83 

5.52 

0.00 

0.25 

3.11 

5.52 

0.00 

0.25 

3.11 

0.20 

2.00 

2.93 

5.46 

0.00 

0.28 

3.42 

5.46 

0.00 

0.28 

3.42 

0.20 

2.00 

3.04 

5.38 

0.00 

0.27 

4.96 

5.38 

0.00 

0.27 

4.96 

0.20 

2.00 

3.08 

5.31 

0.00 

0.29 

5.23 

5.48 

0.00 

0.29 

5.23 

0.20 

2.00 

3.20 

5.03 

0.00 

0.20 

5.48 

5.74 

0.00 

0.20 

5.48 

0.20 

2.00 

3.40 

5.50 

0.00 

0.00 

4.27 

4 . 90 - 

0.00 

0.61 

4.27 

0.20 

2.00 

3.90 

4.88 

0.00 

0.20 

5.71 

4.72 

0.21 

0.78 

5.71 

0.20 

2.00 

4.30 

4.67 

0.00 

0.00 

4.29 

4.50 

0.00 

0.56 

4.29 

0.10 

1.00 

4.44 

4.74 

0.00 

0.00 

4.23 

4.12 

0.00 

0.68 

4.23 

0.10 

1.00 

4.96 

3.58 

0.00 

0.34 

5.74 

5.30 

0.00 

0.11 

5.74 

0.10 

1.00 

5.41 

3.50 

0.00 

0.27 

5.26 

5.38 

0.00 

0.11 

5.26 

0.10 

1.00 

5.87 

4.57 

0.00 

0.05 

6.67 

3.73 

1.01 

0.68 

6.67 

0.10 

1.00 

6.20 

3.26 

0.00 

0.00 

4.22 

4.60 

0.26 

0.00 

4.22 

0.08 

0.50 

6.50 

3.95 

0.00 

0.06 

3.00 

4.83 

0.00 

0.19 

3.00 

0.08 

0.50 

7.20 

3.33 

0.00 

0.00 

6.10 

3.20 

0.21 

0.66 

6.10 

0.08 

0.50 

00 

o 

o 

3.42 

0.00 

0.00 

5.50 

2.56 

0.16 

2.45 

5.50 

0.08 

0.50 

9.50 

3.02 

0.00 

0.83 

6.00 

2.11 

0.00 

2.22 

6.00 

0.08 

0.50 

11.00 

3.64 

0.00 

1.17 

3.01 

1.98 

0.00 

3.78 

3.01 

0.08 

0.50 

13.00 

2.09 

2.23 

0.00 

7.59 

2.09 

2.23 

0.00 

7.59 

0.08 

0.50 
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Table 2. Parameter values for the Saturn proton model (Divine, Ref. 1) 


L a a , a , a b b b b , E E 

2.30 

1.86 

0.00 

0.00 

4.57 

1.86 

0.00 

0.00 

4.57 

1.00 

60.00 

2.32 

2.07 

0.00 

0.00 

4.53 

2.07 

0.00 

0.00 

4.53 

1.00 

60.00 

2.34 

2.17 

0.00 

0.00 

6.09 

2.17 

0.00 

0.00 

6.09 

1.00 

60.00 

2.43 

3.18 

0.05 

0.00 

3.05 

3.18 

0.05 

0.00 

3.05 

1.00 

60.00 

2.52 

2 . 14 

1.85 

0.00 

3.00 

2.14 

1.85 

0.00 

3.00 

1.00 

60.00 

2.67 

3.80 

0.75 

0.00 

3.00 

3.80 

0.75 

0.00 

3.00 

1.00 

60.00 

2.77 

3.37 

1.28 

0.00 

3.00 

3.37 

1.28 

0.00 

3.00 

0.50 

60.00 

2.82 

3.26 

2.63 

0.00 

3.00 

3.26 

2.63 

0.00 

3.00 

0.30 

60.00 

2.86 

3.14 

2.93 

0.00 

3.00 

3.14 

2.93 

0.00 

3.00 

0.30 

60.00 

3 . 10 

0.37 

0.00 

0.00 

3.00 

0.67 

4.87 

6.27 

3.00 

1.00 

60.00 

3 .20 

2.91 

0.29 

0.00 

3.00 

- 0.25 

6.63 

0.00 

3.00 

1.00 

60.00 

3.31 

2.86 

0.02 

0.02 

3.00 

0.29 

8.30 

0.00 

3.00 

0.50 

60 . 00 

3.44 

2.80 

0.98 

0.00 

3.00 

0.46 

7.41 

0.00 

3.00 

0.50 

60.00 

3.71 

2.16 

1.77 

0.00 

3.00 

0.40 

7.10 

0.00 

3.00 

0.50 

60.00 

3.99 

- 0.13 

0.00 

0.26 

3.00 

0.58 

7.70 

0.52 

3.00 

0.50 

60.00 

4 .13 

1.29 

0.00 

0.00 

3.00 

- 0.10 

9.00 

0.15 

3.00 

0.50 

60.00 

4.27 

0.63 

0.46 

0.29 

3.22 

0.16 

8.55 

0.00 

3.22 

0.50 

60.00 

4 .38 

0.61 

5.05 

0.20 

3.12 

- 0.03 

6.93 

0.00 

3 . 12 

0.50 

60.00 

4 .49 

0.47 

5.25 

0.21 

3.55 

- 0.06 

7.15 

0.00 

3 . 55 ^ 

0.50 

60.00 

5.01 

0.68 

1.45 

0.00 

4.41 

- 1.13 

6.73 

6.85 

4.41 

0.50 

0.80 

5.36 

1.77 

0.00 

0.00 

3.00 

- 0.98 

6.78 

8.25 

3.00 

0.50 

0.80 

6.20 

2.32 

3.32 

0.00 

3.08 

0.10 

5.51 

2.99 

3.08 

0.50 

0.80 

6.45 

2.40 

3.34 

0.00 

3.00 

0.06 

5.81 

2.66 

3.00 

0.50 

0.80 

6.70 

2.46 

2.59 

0.52 

4 .72 

0.30 

5.83 

0.12 

4.72 

0.50 

0.80 

7.63 

2.46 

3.55 

0.25 

3.85 

0.81 

4.90 

0.61 

3.85 

0.50 

0.80 

8.40 

2.31 

6.28 

1.73 

3.35 

1.07 

0.48 

1.73 

3.35 

0.30 

0.80 

10.30 

1 . 90 

6.52 

1.81 

4.11 

0.74 

0.00 

0.85 

4.11 

0.30 

0.80 

12.00 

1.24 

0.00 

6.35 

3.94 

1.24 

0.00 

6.35 

3.94 

0.30 

0.80 

13.00 

1.22 

0.00 

6.72 

3.31 

1.22 

0.00 

6.72 

3.31 

0.30 

0.80 
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MAGNETIC EQUATOR 

HIGH LATITUDE 


QUANTITY 

VALUE 

SOURCE 

Radius 

rt*-6.033xl0 7 m 

C084, p. 354 

Gravitation 

GAfj=3.794xlO l# m 3 /s 2 

CE89 

Rotation 

o-8l 0.76 deg/day 

C084, p. 360 

Magnetic field coefficients 

-0.21535 G ^ 
/ 5 -0.01 642 G 
-0.02743 G 
GT--0.0001 G 
G/ --0.00001 G . 

> C084.p. 363. Table II 


Fig. 4. Basic representation of Divine’s Saturn model illustrating the electron and 
proton omnidirectional fluxes and the magnetic field model. 


DIVINE’S CALCULATIONS OF RESIDUALS METHODOLOGY 

A critical step in the Divine model is determining how well it fits the data. 
Unfortunately this is not a straightforward step as the statistical distribution of the data is 
not well defined. Here we review the method used by Divine which assumes an 
approximate log-normal distribution for the fluxes. For each trio of data (observation Rd, 
model R M , background Rb), he computed the following quantities A and Z: 


A = w/n[0,/og(RB/RM)] 
= /og(R D /RM) 

Z = A 2 

= A o 2 [2-(A 0 /A) 2 ] 


for Rd < Rb 
for Rd>Rb 

for |A|<Ao 
for |A|>A 0 


( 4 ) 

( 5 ) 


A 0 was assumed to be 0.3 (giving a factor of ~2 in the ratio of data to model). He then 
computed Y, a modified root-mean-square (RMS) residual, using the weighting factor W: 

Y = [(ZWZ)/(IW)] I/2 (6) 

W is a weighting value, normally assumed to be 1 . The model is iterated until Y achieves 
a minimum value for the assumed database. This process was reviewed by the Board, 
and an alternative will be presented (see below). 
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LA TITUDE AND FLUX VERSUS ENERGY VARIA TIONS 

In addition to the electron and proton omnidirectional fluxes, the model also 
provides the particle pitch-angle variations — that is, how the particle distributions change 
with angle relative to the magnetic field. As will be discussed shortly, the pitch-angle 
distributions assumed by Divine are very simple and may not adequately fit the observed 
distributions. Samples are provided in Figures 5 and 6 in the form of variations with 
magnetic field at the location divided into the magnetic field at the equator for that 
magnetic field line (for a location at the magnetic equator, B = B c , so the ratio is unity; at 
high latitudes, B»B 0 , so the ratio is ~zero). These variations reflect the pitch-angle 
behavior of the model. 

Particle flux variations with energy for the electrons and protons are presented in 
Figure 7 for various L-shell distances (L-shell is approximately the distance in Saturnian 
radii along its equator). As these are integral values, the observation that the flux is 
approximately constant from about 1 to 10 MeV for the protons implies that the model 
assumes little flux between these energies. 




Fig. 5. Latitude variations of the electron omnidirectional integral flux predicted by 
the Divine model in terms of the ratio of the equatorial magnetic field to the observed 
magnetic field. 
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Fig. 6. Latitude variations of the proton omnidirectional integral flux predicted by the 
Divine model in terms of the ratio of the equatorial magnetic field to the observed 
magnetic field. 



Fig. 7. Energy variations of the electron and proton omnidirectional integral fluxes 
predicted by the Divine model for various L-shell values at equatorial and high latitudes. 
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SATRAD MODEL DEVELOPMENT 


STEP # 1 : THE MAGNETIC FIELD AND L-SHELL MODEL 

In order to implement Divine’s algorithms, the B and L coordinates of the 
location must first be provided. As described by R.W. Evans, 1 the determination of the 
(B,L) coordinates can be quite complicated. We have, however, reduced it to a 3 step 
process. First, we define the multipole expansion of the magnetic field and compute the 
field vector and |B| at selected points. This can be quite involved depending on the 
formulation. Table 3 lists the magnetic field components for some of the Saturnian 
models considered. The second step is to trace the field line and integrate the quantity I 
along the field line using the magnetic field model (Ref. 2, Roederer, 1970, Eq. 2.38; see 
also Eq. 4.20): 


1= -B(s)/B m ds 

J s(B m o) V m 


( 7 ) 


The third step requires inverting the I integral to get L corresponding to the point of 
interest. This is given by the following formula from Roederer (Eq. 4.21): 


L 3 R 3 B.. 


= F 


IB.. 


( 8 ) 


«o 7 


Typical results in the form of various L-shell tracings are plotted in Figure 8 for several 
different magnetic field models. It should be noted that in several of the magnetic field 
models, when external field components are included, the magnetic field can go to zero 
or even reverse between 10 to 11 R s , causing the field line tracing program to break 
down. This is a major problem for the radiation model outside ~10 R s . Divine used the 
SZ3 magnetic model (see below) for which L-shells can successfully be traced beyond 
L = 1 0 R s . 


JPL Interoffice Memorandum-5132-04-006, “Planetary Magnetic Field Models for Jupiter, Saturn, 
Uranus, and Neptune; with Field Line Tracing Algorithms (Computer Codes),” January 22, 2004 (JPL 
internal document). 
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Z3 w/o Ring or External Terms 

Z3 w / Ring 

Z3 w / External Terms 


Fig. 8. Magnetic field lines (L-shell contours) for 3 Saturn magnetic field models. 
Field lines are drawn such that they meet at the equator; the lines do NOT represent the 
same L value for all three lines, (see Table 3 for the components for each model). 

The models have been used to estimate the B and L values for the Pioneer and 
Voyager flyby trajectories. These are plotted in B-L space in Figure 9. Besides the 
observation that the field goes to ~0 near 10 R s , it appears from these flybys and from the 
recent Cassini data that the Saturnian magnetic field varies in both UT and Local Time. 
The Saturnian rotation rate has also apparently varied between Pioneer and Voyager and 
Cassini (~6 minutes longer now). Further, the Saturnian field has been difficult to 
represent accurately and has several different representations: Internal, Internal + 
External, or Internal + Ring (we favor Internal + Ring for following field lines) and it is 
(as of this date) not clear which is correct. 

Table 3. Magnetic field coefficients for various models of the Saturn magnetic field. 
(See legend below for abbreviations and acronyms.) 


Coeff 

SZ3 

Z3e 

SP1 

SP2 

SWE 

Z3R 

■ 

Coeff 

Ring 

9i° 

0.21535 

0.21535 

0.21160 

0.21160 

0.21160 

0.21535 

■ 

a 

8.0 

92° 

0.01642 

0.01642 

0.01560 

0.01560 

0.01560 

0.01642 

B 

b 

15.5 

93° 

0.02743 

0.02743 

0.02320 

0.02320 

0.02320 

0.02743 

i 

D_ 

3.0 

Gi° 


- 0.00010 

- 0.000075 

- 0.00012 

- 0.00010 


i 

Hclo 

10 * 

Gi 1 


- 0.00001 



- 0.00001 


1 

1 



Legend 

SZ3: Saturn Zonal Harmonic model 3rd Order 

Z3e: Saturn Zonal Harmonic model with 3rd Order_External Terms 

SP1 : Saturn (Parallel to Earth Sun Line) Voyager model 1 

SP2: Saturn (Perpendicular to Earth Sun Line) Voyager model 2 

SWE: SP1 or 2 with Z3e External Terms 

Z3R: Z3 with Ring 

G m n and g m n are coefficients of zonal spherical harmonic functions. 
The ring model is that of Connerney (Ref. 3). 
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Magnetic Field parameters vs L-Shell 



Fig. 9. Trajectories for the Saturn flybys in magnetic coordinates. Note that Pioneer 
1 1 closely follows the |B cq | contour in to L = 1 .43. 

STEP #2: (B,L) FLUX MODEL IMPLEMENTATION 

The next step in developing the SATRAD model 2 is to implement the model in 
FORTRAN. The process is as follows: 

1 : Look up the coefficients (a n , b n , E 2 , and E 3 ) at discrete L values — linearly interpolate to 
get values at a specific L-shell 

2: Calculate the A„’s as functions of B, B c , and pitch angle a 
3: Calculate I(a,E) and i(a,E) at (B,L) for the A n ’s and E 2 , E 3 


' H. B. Garrett, Jet Propulsion Laboratory Interoffice Memorandum 5130-04-007, “Draft version ot the 
Divine Saturn radiation model code,” April 6, 2004 (JPL internal document). 
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4: Integrate between the pitch angle a’s corresponding to B m and B cr i t (L) at the specified 
position to get the omnidirectional flux. 

The results of the FORTRAN implementation (in color) are overlayed on 
Divine’s original estimates in Figures 10 and 11. Note that there are slight differences for 
the high-latitude components, as we computed the values at the actual B cr j t location rather 
than at the extreme value B e /B=0 assumed by Divine. 
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MAGNETIC SHELL PARAMETER, L (R s ) 


Fig. 10. The FORTRAN modeled omnidirectional fluxes for electrons (solid, colored 
curves) as functions of L for various energies overlaid with Divine's original hand- 
calculated estimates. For the hand-calculations, solid lines correspond to the equator and 
dashed lines to high latitudes-see Fig. 4. 


To test the SATRAD model, it was next used to plot integral flux contours for the 
electrons and protons at 1 MeV and 10 MeV, respectively, in Figure 12. Note the 
complex, multiple belts that appear. 
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MAGNETIC SHELL PARAMETER, L (R s ) 

Fig. 11. The FORTRAN modeled omnidirectional fluxes for protons (solid, colored 
curves) as functions of L for various energies overlaid with Divine's original hand- 
calculations. For the hand-calculations, solid lines correspond to the equator and dashed 
lines to high latitudes-see Fig. 4. 


SATRAD Proton 
10 MeV Integral Flux 
(citF-s )* 1 




West Longitude = 0° 


SATRAD Electron 
1 MeV Integral Flux 
(cnT-s)' 1 

1E6 
1E5 
1E4 
1E3 
1E2 



R s 

Fig. 12. Integral flux contours for the electrons and protons at 1 MeV and 10 MeV, 
respectively, based on the SATRAD computer model (courtesy of I. Jun, JPL). 


STEP #3: COMPARISONS WITH PI I, VI, V2 DATA 

The final step in the SATRAD model development and its review was to compare 
it with the actual data. This is accomplished by computing the residuals between the 
observations and the data. The data sources and instruments and their energy ranges for 
each spacecraft employed in this study are listed in Table 4. The references are listed in 
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the appendices along with additional data references used by Divine. Comparison plots 
for the scanned data and the model predictions are also presented in the appendices (note: 
the data sources used in Divine’s original study are discussed in Ref. 1, which is included 
in the appendices as well). 


Table 4. References, spacecraft, instruments, and energy ranges used in evaluating the 
SATRAD model. 




Electrons 

(MeV) 




Protons 

(MeV) 

REF 

S/C 

Instrument 

Energy 

Range 


REF 

S/C 

Instrument 

Energy 

Range 

VA84 

V02 

LECP 

0.252-0.48 


BC80 

P11 

CPT 

.5-1.8 

PD04 

V02 

LECP 

>0.48] 


VA84 

P11 

GTT 

0.6-3.41 

VA84 

V02 

LECP 

0.48-.853 


CS80 

P11 

CPT-FCD 

>35 

VA84 

P11 

GTT 

>0.56 


SM83 

V02 

CRS 

48-63 

VA84 

V02 

LECP 

0.853-1.2 


KR82 

V02 

LECP 

54-87 

KR82 

V02 

LECP 

>1.5 


SM83 

V02 

CRS 

63-160 

SB03 

V02 

LECP 

>3.24 


VA80b 

P11 

GTT 

>80 

BC80 

P11 

CPT-ECD 

>3.4 


CS80 

P11 

LET 

>80 

BC80 

P11 

CPT-MT 

7-17 


| FM80 

P11 

TRD M3 

>80 

KR82 

V02 

LECP 

10-20 







Acronyms: 

CPT: Charged Particle Telescope (APL) 


-ECD Electron Current Detector 
-MT Main Telescope 
-FCD Fission Cell Detector 
CRS: Cosmic Ray Subsystem 
GTT: Geiger Tube Telescope (U of Chicago) 

LECP: Low Energy Charged Particle Experiment (U of Iowa) 

LET: Low Energy Telescope (U of Chicago) 

TRD M3: Trapped Radiation Detector, channel M3 (UC San Diego) 

RESIDUALS 

To estimate how well the SATRAD model fits the observations, published 
Voyager and Pioneer data were scanned and digitized. Following determination of the 
appropriate geometric factors, the data were next converted to fluxes that could be 
compared to the model predictions. Background values were eliminated from the data 
set. The data and the model comparison plots are presented in the appendices. Based on 
those plots, residuals were estimated using the equations below. 

First assume that the data are distributed log-normally: 


where 


p(r) = 




( 9 ) 


15 


^7 - l°g.o / 

(10) 

obss i 


N 2 


<r=X 

t *- 1 

(11) 


We then define a “standard deviation” or error range: 

STD= 10 ff (12) 

The latter quantity approximates the error range factor we would expect for a true log- 
normal distribution. It is normally used to define the “error range” for Earth’s radiation 
belt models. Figure 13 is a plot of the resulting residuals defined above for the data plots 
considered in this study. The error range is about 2-5.5 (i.e., the typical 
model/observation ratio is inside a range of 2 ±! ) for all the values considered. If the few 
points with residuals outside a range of lO 1 ' (colored blue or dark in Fig. 13) are excluded 
as not being realistic, the overall error range is reduced to a factor of 2.2 (electrons) to 
2.77 (protons) — similar to the factor of 2-3 accepted for Earth trapped-radiation models. 
This method is somewhat different than Divine’s, which limited the residuals above a set 
value of the ratio or below instrument background (see Eqs. 4, 5, and 6). Even so, it 
yields similar results — namely that the SATRAD predictions are typically within 
acceptable limits for a radiation model. 



Fig. 13. The electron and proton residuals calculated for the selected data and the 
SATRAD model are plotted. The total set of data is represented by Blue + Red. The Blue 
points correspond to residual ratios outside the range lO 1 ' . 
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REVIEW BOARD FINDINGS 


The Review Board met on 18 March 2005 at JPL. Table 5 lists the attendees. 
The primary purpose of the board was to determine the suitability of the radiation model 
for evaluating radiation shielding designs for future Saturnian missions. The Board 
determined that SATRAD was currently the best available engineering model for this 
purpose. The secondary purpose was to identify improvements that could enhance its 
value as an engineering tool. Specific questions concerning the model put to the board 
participants are listed in Table 6. To aid in evaluating the Prometheus Saturnian 
Radiation model, reviewers’ responses were requested in the form of “Request for 
Action” sheets. Their responses are presented in Appendix VII. 

Table 5. Saturn Radiation Model Review Board. 


Participant 

In-Person or Call-In 

Institution 

John Cooper 

Call-In 

NASA GSFC 

Joseph Minow 

Call-In 

NASA MSFC 

Giacomo Giampieri 

In-person 

ESA/ JPL 

Tom Armstrong 

In-person 

Fundamental Technology, Inc. 

David Edwards 

In-person 

NASA MSFC 

David Brinza 

In-person 

JPL 

Insoo Jun 

In-person 

JPL 

J. Martin Ratliff 

In-person 

JPL 

Robin W Evans 

In-person 

JPL 

Kent Tobiska 

In-person 

Space Environment Technologies 

Chris Paranicas 

Call-In 

APL 

Charles Barnes 

In-person 

JPL 


Table 6. Specific questions for the Review Board. 


Answer the following questions: 

• Are the high-energy fluxes produced by the model reasonable? 

o Over what energy ranges are Divine’s estimates applicable? 
o Is time-dependence an issue (UT and Local)? 
o What are the problems with the pitch-angle distributions? 

• Are there shortcomings in the magnetic field model? 

o Should external terms be included? 

o Is a ring model more appropriate? 

o How variable is the magnetic field (UT and Local)? 

• Does the particle model adequately represent the Pioneer and Voyager data? 

o At the end of the day, does it fit the data? 

o Is the method used to evaluate the ability of the model to fit the spacecraft 
data reasonable? 

• Is the model a useful engineering design tool as formulated? 
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MODELING ISSUES 

Although the review brought to light concerns with the original model and its 
current implementation. Divine was well aware of various problems in fitting the original 
Pioneer and Voyager data and in the fitting process itself. These can be listed as follows. 

• Actual data points versus lines (or fits) 

• Inbound/Outbound observations mismatch 

• Raw versus corrected count rates (i.e., dead-time corrections) 

• Omni versus spin-modulated responses 

• Species confusion and alternation 

• Detector response functions 

• Time-averaged versus spin-averaged count rates 

• Background count-rate levels 

• Mistakes in computing the differential spectra (e g., AI/AE replacing i) 

• Logarithmic versus linear scales 

• L-shell (a magnetic field coordinate) versus R (radial distance) scales 

• Left- versus right-hand scales (often not clear which goes with which) 

• Factors used to spread the profiles 

• Pitch-angle versus latitude distributions (see Figures 5 and 6 for example) 

In an attempt to address these issues, the original model developed by Divine was 
reviewed immediately after its initial issuance in March 1990 (almost exactly 15 years 
ago) by outside reviewers. As their comments (based on written correspondence to the 
Project) closely follow those at the March meeting, they are summarized here. 

Richard Selesnick of the California Institute of Technology, for example, noted 
that the pitch-angle data he was aware of changes its dependence to isotropic near 4 R s — 
this implies the model is too anisotropic outside this distance. He suggested that there 
was a need to model each encounter separately to test if there are any significant time 
variations. 

S. M. Krimigis of the Applied Physics Laboratory (APL), who was principal 
investigator for the Low-Energy Charge Particle (LECP) [experiment] instrument (a main 
source of data for the model), also reviewed the data. He indicated that the pitch-angle 
dependences were not in agreement with the APL observations of protons at L~2.67 and 
electrons at L-4.5-5. He believed that the “protons only” assumption that Divine made in 
interpreting ion flux data was dubious at energies <1 MeV (e g., oxygen ions may be 
present). Although Divine ignored the LECP ion data at various locations, Krimigis 
believed that the LECP ion spectra at L = 7. 15 R s in the inner magnetosphere were 
reliable and should be included in the model. He also indicated that he felt that two 
instrument detectors (the S 180B and MD 80 detectors) used by Divine were 
contaminated by electrons and may have had high proton fluxes which made these data 
unreliable. 
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The most extensive review was carried out by John Cooper, one of the current 
reviewers (see Table 5). He wrote that there were background contamination concerns in 
several data sets. Sub-MeV data should only be used with caution inwards of the Tethys 
orbit (4.88 R s ) as penetrating trapped radiation in the inner magnetosphere mimics the 
response to sub-MeV particles. Electrons and the 100 MeV proton-induced background 
overwhelms the 0.5-10 MeV proton data in the Saturnian proton belt. The Voyager 2 
LECP responded to background from higher-energy electrons while the radial profiles for 
sub-MeV electrons were similarly affected by high-energy background. Only the 100 
Mev proton curves were good inside 5 R s . Lower energies were contaminated and 
represented only an upper limit. Cooper also had various concerns with the pitch-angle 
distributions — he has since suggested that the model needs to use the more flexible pitch- 
angle distributions of the Divine-Garrett Jupiter radiation model that treat the loss cone 
better. It should be possible to include more Pioneer 1 1 (unpublished) data to determine 
the electron anisotropies. Likewise, the model needs to improve its description of the 
proton anisotropies. In fact, the model requires B c = 0 for A n = b n and allows only 
isotropic fluxes at high latitudes — it needs to be updated with Cassini data. 

The pitch angles are more complex than isotropic inside ~4 R s (probably not 
important for radiation modeling). The 63-160 MeV proton anisotropies have strong 
radial variations so that the uniform isotropy assumed by the model in the inner region is 
incorrect (as a specific concern, what is source of 100 MeV anisotropy profile at 7.63 
R s 7). There is probably real-time variability in the MeV electrons that the static Divine 
model does not account for. There are a few specific glitches in the model — namely the 
10 MeV peak at 2.5 R s and the notch at 2.65 R s are in error (see Figure 4) — need to 
review the Northrop and Fillius Pioneer 1 1 data. Finally, the model needs to show that 
electrons dominate in satellite regions — protons elsewhere. The 0-20 MeV electrons 
determined from the profile in Krimigis et al., 1982 (see Ref. KR82 in appendices) may 
be an order of magnitude too high at 2.7 R s in the model. 

REVIEW BOARD COMMENTS 

In addition to the above, the Review Board comments on the SATRAD model can 
be roughly divided into 4 groups: 

1) Pioneer and Voyager particle data (Note: these data are still of value as they cover 
energy ranges and locations not covered by Cassini): 

a) Need to re-evaluate using data now in the Planetary Data System (PDS). 

b) Need to re-evaluate detector background sources and calibrations in the 
light of current understanding of the instruments. 

c) Potentially important local-time (relative to Sun) and real-time “weather” 
variations — difficult to separate! 

d) Species identification requires more careful consideration — ions versus 
electrons, protons versus oxygen ions; etc. 

e) Some spurious features in model (see J. Cooper comments). 

f) Pitch-angle distributions in particular need to be reconsidered and a more 
flexible fitting scheme employed. 


19 


2) Cassini data: 

a) Need to incorporate latest Cassini data. Unfortunately, only data from the 
Saturn Orbit Insertion (SOI) orbit may be useful for critical pitch-angle 
information as instrument may now be “stuck” in one position 

b) Cassini found an inner belt (see Figure 14) in a region that was previously 
unexplored, and thus not currently included in the model. 

3) Magnetic field: 

a) Uncertainties in magnetic field make field line tracing difficult at ~10 R s 
as field approaches zero there. 

b) As for particles, potentially important local time (relative to Sun) and real 
time (“weather”) variations that are difficult to separate! 

c) Cassini found 6 minute slowdown in Saturn rotation and a highly variable 
magnetic field. 

d) Need to include effects of magnetopause and ring current in magnetic field 
model in outer magnetosphere (~10 R s ). 

4) Overall model package 

a) Need to standardize and have peer-reviewed in literature. 

b) Need to make code/engineering model available for Cassini and other 
programs outside JPL. 


REVIEW BOARD RECOMMENDATIONS 


Based on the preceding, the Review Board recommendations as to areas that 
should be considered in both using the model and in updating it in the future can be 
concisely summarized as follows: 

• Calibration/background issues with the original Pioneer and Voyager data 

need to be taken into account. 

• The pitch-angle distributions need to be reconsidered. 

• There are concerns over the apparent variability of the Saturnian 

magnetosphere with respects to local solar time and absolute time. 

• Cassini data should be actively pursued to improve the SATRAD magnetic 

field model and flux estimates. 

• The model should be standardized, made available to the general community, 

and formally reviewed in the literature. 


20 



Fig. 14. Energetic neutral atom (ENA) image obtained during Cassini’s passage over the 
rings showing an inner radiation belt not included in the current model (A) Ion and 
Neutral Camera (INCA) image for 20-to 50-keV/nucleon ENAs. Bright region above 
Saturn’s limb is caused by ENAs produced by charge-exchange collisions between main 
radiation belt ions and near-equatorial gas. The band of emission above equator on disk is 
produced by the same ENAs from the main belt being stripped in Saturn’s exosphere 
between the inner edge of the D ring and the cloud tops, trapped there temporarily as 
energetic ions and then reemitted as ENAs that INCA can image. The image is bounded 
by the limits of the INCA field of view (dotted arcs). (B) Schematic of charge 
exchange/stripping process that begins as ENA emission from the main belt and produces 
ENA emission from Saturn’s exosphere. (Reprinted with permission from Krimigis et al., 
“Dynamics of Saturn’s Magnetosphere from MIMI during Cassini’s Orbital Insertion,” 
Science, Vol. 307, pp. 1270-1273, February 25, 2005. Copyright 2005 AAAS). 

CONCLUSIONS 

In the early 1990’s, Neil Divine developed a Saturn proton/electron radiation 
model based on data from Pioneer 11, Voyager 1, and Voyager 2. That model has been 
turned into a FORTRAN computer code for analyzing the radiation environment at 
Saturn. The code, SATRAD, uses Divine’s formulation to compute the electron and 
proton fluxes at a specified B and L coordinate. A separate model of the Saturnian 
magnetic field is used to compute the B and L coordinates along a spacecraft trajectory. 
While originally intended as an engineering tool to model radiation doses to Cassini 
during its design, the current model implementation was used to estimate the doses for 
Prometheus missions to Saturn and will be used to compare with the Cassini 
measurements when they become available. To validate the SATRAD model, its flux 
predictions have been compared by the authors with published Pioneer and Voyager data. 
Those comparisons allow an estimate of the “accuracy” of the SATRAD model by 
providing an error range. Within the limits of the current study, the original estimate of 
Divine has been validated to a factor of 2-3 and is within current expectations for 
radiation modeling. While agreeing to the general usefulness and validity of SATRAD 
for radiation modeling, the Review Board found several areas that should be considered 
in both using the model and in updating it. These ranged from calibration/background 
issues with the original Pioneer and Voyager data to very real concerns over the apparent 
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variability of the Saturnian magnetosphere. Finally, Cassini data, when they become 
available, should allow improvements in the SATRAD magnetic field model and flux 
estimates and should be actively pursued. The reviewers requested copies of the current 
model for comparisons with the Cassini data and recommended that the model be 
standardized, made more readily available, and formally reviewed in the literature. This 
document was written to address those issues. 

TEXT REFERENCE 

1. Divine, Neil, JPL Interoffice Memorandum-5217-90-029, Numerical Models for 
Electron and Proton Distributions in Saturn's Radiation Belts, Jet Propulsion Laboratory, 
Pasadena, California, 28 February 1990. 

2. Roederer, J. G., Dynamics of Geodynamically Trapped Radiation, Springer- Verlag, 
Berlin, 1970. 

3. Connemey, J. E. P., M. H. Acuna, and N. F. Ness, “Currents in Saturn’s 
magnetosphere”, J. Geophys. Res., V88, A1 1, pp. 8779-8789, 1983. 

4. Krimigis, S. M., D. G. Mitchell, D. C. Hamilton, N. Krupp, S. Livi, E. C. Roelof, 
J. Dandouras, T. P. Armstrong, B. H. Mauk, C. Paranicas, P. C. Brandt, S. Bolton, 
A. F. Cheng, T. Choo, G. Gloeckler, J. Hayes, K. C. Hsieh, W.-H. Ip, S. Jaskulek, 
E. P. Keath, E. Kirsch, M. Kusterer, A. Lagg, L. J. Lanzerotti, D. LaVallee, J. Manweiler, 
R. W. McEntire, W. Rasmuss, J. Saur, F. S. Turner, D. J. Williams, J. Woch, “Dynamics 
of Saturn’s Magnetosphere from MIMI During Cassini’s Orbital Insertion”, Science , Vol. 
307, pp. 1270-1273, February 25, 2005. 


22 


APPENDICES 


23 


APPENDIX I. Table Al. Electron Model Data Sources ( Divine , Ref. 1) 
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NOTES TO TABLE A1 


[a] The profiles appear in order of decreasing threshold energy E a (see column 5 above and note [c], 
below); a profile lacking an entry in column 1 is used in the model only for trajectory points at 
which no value for count rate R (see note [d]) is available from the profile identified on the line 
immediately above. 

[b] Second entry in this column is the curve placement, counting down from the top of the figure 
cited. 

[c] Where two entries appear, the response (next column) is assumed constant between them, and zero 
elsewhere; a single entry represents the smaller (threshold E a ) of these two energies, and the other 
(E b ) is assumed infinite. 

[d] Units cm 2 assumes R is a count rate (units s' 1 ); exceptions are cited in the notes, and the entry in 
this column is always the ratio of the profile quantity R (count rate, intensity, flux, etc.) to the 
local omni flux J (see text following Eq. 7). Parentheses enclose powers of 10. 

[e] This entry is estimated from the noise or limit of the profile in the figure, in the same units as R in 
the preceding column (see note [d]). 

[f] Energies from profile label; conversion is sG in cm 2 sr (KR83, p. 8874, Table 1) divided by 47i sr. 

[g] Energy from figure label, profile is flax J, no conversion needed. 

[h] Energy from source (FM80, p. 5805, Table 1), profile is current; conversion is as discussed in 
Appendix B (last paragraph). 

[i] Energies from profile label, conversion is 0.43 cm 2 sr (SI80b, p. 5757, Table Al) divided by 47i sr. 

[j] Energy from profile label, conversion is sG=0.12 cm 2 sr (KR83, p. 8874, Table 1) divided by 4n sr 
and by an additional factor 10 to secure agreement with a corresponding profile (KA82, p. 1 144, 
Fig. 2, curve 3 from top) and with Pioneer 1 1 detectors at nearby energies. 

[k] Energies from figure caption (or label), profile is differential intensity i. conversion is (47i sr)" 1 (Eb- 

Ea)' 1 . 

[l] Same as note [f], but conversion also confirmed by comparison with spectrum at L = 9.56, as 
discussed in Appendix B (paragraph 4). 

[m] Energies from source (VA80a, p. 416, Table 1), conversion is Q from same table. 

[n] Energies as in note [m], profile is integral intensity I, conversion is (47t sr)" 1 . 

[o] Energy from profile label, conversion is (R/J) 3 (Ri/R 3 ), where subscript 1 is for channel El, 
subscript 3 is for channel E3, Ri and R 3 are the peak count rates at L = 2.7, and conversion (R/J) 3 
is derived in note [p]. 

[p] Energy from source (FM80, p. 5805, Table 1), conversion is 0.0057 cnrsr (same table) divided by 
4ir sr. 

[q] Energies from source (F180, p. 426, Table 1), profile is differential intensity i, conversion is 
(47t sr )"‘(E b — Ea)" 1 . 
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APPENDIX II. Table A2. Proton Model Data Sources (Divine, Ref. 1 ) 
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NOTES TO TABLE A2 


[a] The profiles appear in order of decreasing threshold energy E a (see column 5 above, note [c] 
below); a profile lacking an entry in column 1 is used in the model only for trajectory points at 
which no value for count rate R (note [d]) is available from the profile immediately above. 

fb] Second entry in this column is the curve placement, counting down from the top of the figure 
cited. 

[c] Where two entries appear, the response (next column) is assumed constant between them, and zero 
elsewhere; a single entry represents the smaller (threshold E a ) of these two energies, and the other 
(E b ) is assumed infinite. 

[d] Units cm' 2 assumes R is a count rate (units s' 1 ); exceptions are cited in the notes, and the entiy in 
this column is always the ratio of the profile quantity R (count rate, intensity, flux, etc.) to the 
local omni flux J (see text following Eq. 7). Parentheses enclose powers of 10. 

[e] This entry is estimated from the noise or limit of the profile in the figure, in the same units as R in 
the preceding column (see note [d]). 

[f] Energy from profile label, conversion is geometric factor in cm 2 (FI80, p. 426, Table 1). 

[g] Energy from source (FM80, p. 5805, Table 1), conversion is geometric factor in cm 2 (same table; 
note disagreement with FI80, p. 426, Table 1). 

[h] Energy from source (VA80a, p. 416, Table 1), conversion is Q from same table. 

[i] Energy from source (VA80b, p. 5682, Table 2), conversion is Q from same table. 

[j] Energy as in note [h], conversion from source (VA83, p. 6912, Eq. 1). 

[k] Energies from figure caption (or label), profile is differential intensity i, conversion is (47~ sr) 1 

(Eb-Ea)-'. 

[l] Energies from profile label, conversion is sG in cm 2 sr (KR83, p. 8874, Table 1) divided by 4 n sr. 

[m] Energies as for preceding profile, conversion is (4n sr)' 1 (E b -E a ) _1 (R3A3) where logR 3 =3.75 and 
log i 3 =1.21 are peak rates from this and the preceding profile. 

[n] Energy from profile label, conversion from source (SI80b, p. 5757, Table Al), weight is 2.0 to 
compensate for indistinguishability of the inbound and outbound profiles. 

[o] Energies from figure caption, profile is integral intensity I, conversion is (4 n sr) 1 . 

[p] Same as note [k], but conversion also confirmed by comparison with spectrum at 7.81 R*. as 
discussed in Appendix B (paragraph 4). 

[q] Energies from profile label, conversion is geometric factor in cm 2 sr (SI80b, p. 5757, Table Al) 
divided by 4n sr. 
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APPENDIX IV. Data Comparisons between SATRAD and Pll, VI, AND V2 


Fig. A1 (next page and following) SATRAD model predictions versus PI 1, VI, and 
V2 data. Specific frame, reference (see Appendix III), and figure number for each 
reference are as follows: 
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V2e0.48<E<0.853 MeV 



Pll e E>0.56 MeV 



V 2 e 0.853<E<1.2 MeV 



V2 e E>1.5 MeV 
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Electron Flux (cm-2 s-1) 


V 2 e E>3.24 MeV 



Pll e 7-17MeV 


L-Shell 


P11 aF >3.4 MeV 




x 


8 


■ Pll out (cm-2 
s-1) 

■ Pll E >3.4 
model 



V2 e 10<E<20 MeV In 



V2 e 10<E<20 MeV Out 


l.E+03 


l.E-02 


l.E^OO ■ 


Pll p 0.5<E< 1.8 MeV 


J 









(cm-2 s-1) 


X 



Model 10<E<20 


to 











1 

h 




11 

V 




2 

4 

L-S 

6 

.hell 

8 ^ 

10 



L-Shell 
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Pll p 0.6<E<3.41 MeV in 



L-S^ell 


Pll p E > 35 MeV 


E 

| 

2 1 E+02 
o 

£ 


s 

i 




L 






(cm 

i-2 s-1) 

5 Mev 


1 

r 


■■^■Pll > 3 










1 

5 





6 

L Shell 


l.E-02 
l.E-01 
^ 1.E-00 

M 

Ft 

g l.E-01 
u 

9 


l.E-02 


l.E-03 


l.E-04 


l.E-05 


1 






4 



^■^2 (cm-2 s-1 sr-1 MeV- 
1) 

Model p 48 < E < 63 

J 

ft 





i 


(L 





1 1 

\ 




1 



i; 

7 





10 


Pll p 0.6<E<3.41 MeV Out 


O 1.E+-00 



■ Pll (cm-2 s-1) 

■ Model p 0.6 < E < 3.41 


14 


6 

L- Shell 


V2 p 48<E<63 MeV in 



a.E-i-oi 


d -E-POO 


l.E-03 


1 





l 





1 

1 



■■BHV2 (cm- 
Model p 

2 s-l)-ln 
54<E<87 



j 




i 












1 



1 

1 

» 



1 


L-sftell 
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l.E+05 

l.E+04 

l.E+03 

\.E+02 

<N 

i 

■3.E+01 

J2 

U. 

Bl.e+oo 

I 

l.E-01 

l.E-02 

l.E-03 


V2 p 54<E<I 

B7 MeV 





"1 



HV2 

ou 

"■Mr 

. (cm-2s-l)- 
t 

iflel n 

n 

A 










A 

k 




1 

F 1 

1 






1 



1 

9 

1 

L 




L-#ie!l 


10 


V2 d 63<E<160 MeV Out 



l.E+02 


V2 p 63<E<160 MeV in 


V2 (cm-2 s-1) 

Model p 63 < E < 160 


l.E-01 


l.E-02 


l.E-03 


l.E-04 



Pll p E>80 MeV 


— 

A 




O' 

L 


l^^HPll (cm-2 s-1) 

Model E>80 _ 


■ 





n 













22 





L-^iell 


10 


Pll p E>80 MeV M3 

1 t *05 



lE^OO — 

2 4 & « 1 ° 

L-Shell 


Pll P E> 80 MeV Lin2 LET 



L-Shtll 
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APPENDIX V. SA TRAD Listings 


Saturn Radiation Program 

The following is a FORTRAN listing of the SATRAD model. It consists of a 
main calling program, “Smainrxx,” which reads in the trajectory and magnetic field data 
along the trajectory (see “traj.in”) and prepares various output files. It also establishes 
the energy channels to be analyzed. This program can be easily modified to handle a 
variety of input formats. The second primary component is subroutine SATRADxx and 
its subsidiary subroutines. This component should be kept as is. 


Program Smainrxx 

c Smainrxx. f is an example of a main calling program for 

c SATRAD. f SATRAD. f calculates the Saturnian electron and proton fluxes, 

c 

c The SATRAD Program is documented in JPL Publication 05-9, October 2005, 

c "Saturn Radiation (Satrad) Model", H.B. Garrett, J.M. Ratliff, 

c and R.W. Evans Algorithms are those of Neil Divine, as documented in 

c N. Divine, "Numerical Models for Electron and Proton Distributions 

c in Saturn's Radiation Belts", JPL IOM-5217-90-029, 28 February, 1990. 

c Final coding by Garrett and Evans, based on subroutines created by Ratliff 

c ca. 1992, with I/O changes and some restructuring courtesy of Janet 

c Luhmann, UC Berkeley, ca . 1994. 

c H. Garrett 


c Smainr is an example of how to provide the set up variables 

c necessary for the Saturnian radiation model in the "Satrad" package. The 
c inputs and common variables are listed below. Note that the primary 
c input is observer position, species (electron or proton) , and energy, 
c The Satrad package returns the electron and proton differential 
c and integral omindirectional energy flux. 

c Variables : 

c [variable names are enclosed in single quotes... 'variable'.] 


c 

c 

c 

c 

c 

c 


INPUTS : 

'EE(n) ' = Vector of energies at which fluxes are to be 

estimated; "n" is the number of energies to be computed. 
Units are MeV. 

' E ' = Input energy at which fluxes are to be 

estimated by Satrad. Units are MeV. 


c ' N ' 

c 


= Number of integration steps in Eq. 10, Ref. 1 

(typically, N = 7) 


Coordinate System: 
'ALAT' 


'RS' 


' WLONG ' 


Input latitude at which fluxes are to be 
estimated by Satrad. Units are degrees 
(north is positive, south negative) . 

= Input radial distance at which fluxes are to be 
estimated by Satrad. Units are Saturnian radii 
(1 Rj = 60,330 km) . 

= Input longitude at which fluxes are to be 
estimated by Satrad. Units are degrees 
(longitude is positive westward) . 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OUTPUTS : 
' OFLX ( 4 ) 


= Output vector from Satrad (Note: all values are 
log base 10) : 


OFLX(l) =J(CM**2-S) **-l 

OFLX (2) =DJ/DE ( CM* *2 -S -MEV) **-l 

OFLX (3) =J (CM**2-S) **-l 

OFLX (4 ) =DJ/DE (CM**2 -S-MEV) **-l 


for electrons 
for electrons 
for protons 
for protons 


where : 


J = omnidirectional integral flux as given by 
Eq. 4, Ref. 1 
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c DJ/DE = omnidirectional differential flux 

c 'OOF(4,n)' = Output matrix of (10**OFLX) values at energies l,...,n. 

c 'VX' = Output vector such that: 


c 

V(1) 

= 

Bx component of magnetic field (units of G) 

c 

V (2) 

= 

By component of magnetic field (units of G) 

c 

V (3) 

= 

Bz component of magnetic field (units of G) 

c 

V (4) 

= 

B magnitude of magnetic field (units of G) 

c 

V (5) 

= 

L value of field line through point 

c 

V (6) 

= 

Cutoff latitude of field line through point 

c 



(cutoff means top of atmosphere) 

c 

V(7) 

= 

Be, amplitude of field line at equator 

c 

V (8) 

= 

Be, cutoff field strength along field line 

c 

V (9) 

= 

B, amplitude at mirror point for given pitch 

c 



angle 

c 

V(10) 

= 

theta, magnetic co- latitude (in degrees) of 

c 



point 

c 

V(ll) 

= 

OFLX(l) 

c 

V ( 12 ) 

= 

OFLX (2) 

c 

V (13) 

= 

OFLX (3) 

c 

V ( 14 ) 

= 

OFLX (4) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


' ALE ( 21 ) ' 

1 AAE (14,21) 1 
' ALP (28 ) ' 

' AAP (14,28) 1 
1 BB (14 , 19) 1 


BB ( 1 , 1-19) 
BB (2,1-19) 
BB (3, 1-19) 


= L shell values for electron model (Table 5, Ref. 1) . 

= Parameter values for electron model (Table 5, Ref. 1) . 

= L shell values for proton model (Table 6, Ref. 1) . 

= Parameter values for proton model (Table 6, Ref. 1) . 

= Matrix of critical latitude and B field and of B at 

the magnetic equator such that (note: BB(4-14 # 19) are 
not used 11): 

= Equatorial field strength at L values 
given in vector ALO (see Satrad) 

= Cufoff field strength at L values 

given in vector ALO (see Satrad) 

= Cufoff latitude at L values given in vector 
ALO (see Satrad) 


DIMENSION EE (13) , OFLX (4) # OOF (4 ,13) ,OOFT(4,13) ,VX(14) 
character *80 aa 
COMMON/ RADAT/VX 

c Input array values for energy steps (EE) 

DATA EE/.l, . 2 , .3, .5, 1.0, 2. 0,3. 0,5.0, 10. 0,30., 100. , 300 . , 1000.0/ 

c Default values for west longitude (WLONG = 0; the magnetic and 

c spin axes are assumed to be aligned) and the number of 

c angular integration steps (can be as high as ~98) . 

DATA WLONG, N/0 . 0 , 14 / 

c Main output data files 

c For this exercise, we will use TIME = 1 minute = 60 s 

c e . g . , time = 60 . 

c Here we will read time from input file. . . 

DO 10 J=1 , 4 
DO 10 IE=1 , 13 
10 OOFT ( J, IE) = 0 

OPEN ( 2 , FILE= ' SOrbitf luence . dat ' ) 
write (2,*) ’Divine, IOM 5217-90-029, 02/28/90’ 
write (2 , * ) 1 ' 

write (2 , * ) ’ FLUENCE SPECTRA' 

OPEN ( 5 , FILE= ' Sf luence . dat ’ ) 

write (5,*) 'Divine IOM 5217-90-029, 02/28/90' 
write (5 , * ) ' ' 

write (5,*) 'FLUENCE SPECTRA' 
write (5, *) 1 Electrons 

x Protons ' 
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Integral 
(cm** -2 ) 


Differential 


Integ 

(cm** 


write (5,*) 'Energy 
xral Differential' 
write (5,*) ' (MeV) 
x-2) (cm**2 MeV)**-l' 


(cm* *2 MeV)**-l 


OPEN (4 , FILE= ' Sf luxint . dat ' ) 

write (4 ,*) 'Divine IOM 5217-90-029, 02/28/90' 
write (4,*) 'Energy (MeV) levels E(l) to E(13) are-' 
write (4 , 103 ) (EE (i ) , i=l,13) 
write (4,*) ' ' 

write (4 , * ) ' INTEGRAL FLUX SPECTRA' 

write (4,*) Each line contains the following sequence of informatio 


write (4 , 104 ) 

write (4 ,*) 'Time (JD) Range(Rs) Lat(deg) WLong(deg) B(gauss) Beq ( 
L-shell Blat(deg) Se(E(l) Se(E(2)) Se(E(3)) Se(E(4)) 

X ?f\ E i 5 li/-. S ?» E(6) ) Se(E<7,) Se (E ( 8) ) Se (E ( 9) ) Se(E(10)) Se(E(l 

xl ) ) Se(E(12)) Se(E(13)) Sp(E(l) Sp(E(2)) Sp(E(3)) Sd(E(4)) Sd 

x) (E Sp(E(12n E Sp(E(13?r (7) } SP(E(8>) S P (E<9)) Sp(E(10)) Sp(E(ll) 


write (4 , * ) 1 ' 


OPEN ( 3 , FILE= ' Sf luxdi f f . dat ’ ) 

write (3,*) ’Divine IOM 5217-90-029, 02/28/90' 

write (3 ,*)' Energy (MeV) levels E(l) to E(13) are* ' 

write (3 , 103 ) (EE (i ) , i = l,13) 

write (3,*) ' 1 

write (3,*) 'DIFFERENTIAL FLUX SPECTRA' 

write (3 ,*)' Each line contains the following sequence of informatio 


write (3 , 105) 

write (3 ,*)' Time (JD) Range(Rs) Lat(deg) WLong(deg) B(gauss) Beq ( 
, 1 L ' she11 Blat(deg) Se(E(l) Se(E(2)) Se(E(3)) Se(E(4)) 

xSe (E (5) ) Se (E (6) ) Se(E(7)) Se(E(8)) Se(E(9)) Se(E(10)) Se(E(l 

WP/cf? (E<18)) S P( E <1) Sp (E (2) ) Sp ( E ( 3 ) ) Sp (E ( 4 ) ) Sp 

X) Sp(E(12n Sp(E( 13 n E ?>) SP(E<8)> S P (E(9)) Sp (E (10 ) ) Sp(E(ll) 


write (3, *) ' 1 


Loop to setup position (rs , lat , wlong) 

Read from file traj .in. Format is julian date, radial distance 
m Rs, latitude (deg), and west longitude (deg) 


open (8 , f ile= ' traj .in') 

Header files 

read (8, 101) aa 
write (*,*) aa 
read (8, 101 ) aa 
write (*,*)aa 
read (8, 101) aa 
write (*, *)aa 

Lines of data 

read (8 , * ) NN, time 
write (* , * ) NN, time 
read (8, 101 ) aa 
write (*, *) aa 

do 2 irj =1, 1000 


read (8 , *) date, rs, alat , wlong, al , B, Be, Be 

Check for end of orbital segment 
if (date .eq. 1.0) go to 22 
Check for end of data 

if (date .eq. 2.0) go to 3 


Loop to setup energy channel 


DO 1 IE=1 , 13 
E=EE (IE) 

c Main call to Satrad library 

c CALL SATRAD ( E , ALAT # Rs , WLONG f N , OFLX ) 

Cal 1 Satradxx (E # B , Be , Be , AL # ALAT f Rs f WLNG # N # OFLX ) 

c Check loop to catch "-999" zero value flags (actually the 
c -999 represents 0.0 value as all OFLX values are log base 10 and 

c l0**-999=0 is valid...). Note that OFLX and OOF are log base 10. 

DO 1 J=1 , 4 

IF (OFLX (J) .EQ.-999.) OOF ( J, IE) =0 . 0 
IF (OFLX ( J) . EQ . - 999 . ) GO TO 1 
OOF ( J, IE) =10 . **OFLX ( J) 

OOFT ( J, IE) =OOFT ( J # IE) +OOF ( J, IE) *TIME 

1 CONTINUE 

c Output : 

b=vx (4 ) 
be=vx (7) 

AL=vx ( 5 ) 

BLAT=90 . -VX (10) 

WRITE (4 # 100) date, rs , alat , wlong, b, be, al , blat , (oof (1 , ie) , 
x ie=l, 13) , (oof (3 , ie) , ie=l, 13) 

WRITE (3,100) date , rs , alat , wlong, b, be, al , blat , (oof (2,ie) , 
x ie=l, 13) , (oof (4,ie) ,ie=l, 13) 
dater=date 
go to 2 

22 CONTINUE 

write (2 , * ) dater 
write(2,*) ! Electrons' 

write (2 , * ) 1 Energy Integral ' 

write (2,*) ' (MeV) (cm**-2) ' 

do 5 ie = 1,13 

5 write (2 , 102 ) ee (ie) , ooft (1 , ie) 

write (2 , * ) ' Protons ' 

write (2 , * ) 1 Energy Integral ' 

write (2,*) ' (MeV) (cm**-2) ' 

do 6 ie = 1,13 

6 write (2 , 102) ee (ie) , ooft (3 , ie) 
write (2/ *) * ’ 

2 continue 

3 CONTINUE 

do 4 ie = 1,13 

4 write (5 , 102 ) ee (ie) , (ooft (j,ie) , j = l,4) 

100 format (f 10. 4, 7 (lx, e9. 3) , 13 (lx, e9 . 3) , 13 (lx, e9 . 3) ) 

101 format (la80) 

102 format ( f 9 . 2 , 4x, e9 . 3 , 3x, e9 . 3 , 15x, e9 . 3 , 4x, e9 . 3 ) 

103 format (13 (flO .2) ) 

104 format (14 Ox, ' Electron Flux (cm x 2- s) A -1 ' , 120x, 
x' Proton_Flux (cm x 2-s) : * : -l) ' ) 

105 format (140x, ' Electron_Flux (cm A 2- s-MeV) A - 1 ' , 115x, 
x’Proton_Flux (cm A 2-s-MeV) A -l) ' ) 

CLOSE ( 2 , STATUS= ' KEEP ' ) 

CLOSE ( 3 , STATUS= ' KEEP ’ ) 

CLOSE ( 4 , STATUS= ' KEEP ' ) 

CLOSE ( 5 , STATU S= ’ KEEP 1 ) 

pause 

STOP 

END 


c 'VX' = Output vector such that: 
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on OOOOOOOOOQOOO 


V ( 1 ) = Brad radial component of magnetic field (units of G) 

V (2 ) = Btheta component of magnetic field (units of G) 

V (3 ) = Bphi component of magnetic field (units of G) 

V (4 ) = B magnitude of magnetic field (units of G) 

V ( 5 ) = L value of field line through point 

V ( 6 ) = Cutoff latitude of field line through point 

(cutoff means top of atmosphere) 

V ( 7 ) = Be, amplitude of field line at equator 

V ( 8 ) = Be, cutoff field strength along field line 

V (9 ) = B, amplitude at mirror point for given pitch 

angle 

V (10) = theta, magnetic co-latitude (in degrees) of point 
V(ll) = OFLX(l) 

V ( 12 ) = OFLX ( 2 ) 

V ( 13 ) = OFLX ( 3 ) 

V ( 14 ) = OFLX ( 4 ) 

************************* 

Contains RDD/DMU corrections. 

SUBROUTINE SATRADxx (E , B , Be , Beq, AL , ALAT, RS , WLONG , N , OFLX) 

DIMENSION VX ( 14 ) ,XX(100) ,CUR(4,100) ,CUR1(100) ,OFLX(4) 

COMMON/ RADAT/VX 
C** Corrected 09/2000 (more accurate) HBG 
RDD=. 017453292 

ALPHA=0 
Vx (1 ) =0 
Vx (2 ) =0 
Vx (3 ) =0 
Vx (6 ) =0 
Vx (4 ) =B 
Vx (5 ) =AL 
Vx(7) =Beq 
VX (8 ) =BC 

BM=B/ (SIN (ALPHA*RDD) **2) 

Vx(9)=Bm 

Theta=90 . -alat 
Vx (10 ) = Theta 

c If the L shell is less than 2.3 or greater than 13, then 
c the model gives 0 values. As the data are log base 10, a 
c common value of -999 is assumed to indicate "0" . 

IF (AL . LT. 2 . 3 .or. AL.GT.13.) GO TO 55 

c If the magnetic field amplitude B is greater than BC, then 
c the particle is in the loss cone and its flux is set to "0". 

IF (B . GE . BC) GO TO 55 

c ’ dMu' = differential step in pitch angle used in the integral, 
c It is given by dividing the critical pitch angle by N-l intervals. 

BMUC=SQRT (1 . -B/BC) 

DMU=BMUC/ (N-l) 

c** Corrected for numerical inaccuracy 09/2000 HBG 
dmu=dmu- . 000001*dmu 

c This is the main loop to set up the pitch angle integrations. 

DO 1 IA=1,N 

XX ( I A) = FLOAT ( I A- 1 ) *DMU 
ALPHA=ACOS ( FLOAT ( IA- 1 ) *DMU) /RDD 

c ’ALPHA' = pitch angle of particle. 

CALL SRADBL ( B , B c , Beq , AL , ALAT , RS , WLONG , ALPHA , E ) 


c 

c 

c 
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o n o n o 


CUR (l.IA) =VX(11) 

CUR (2 , I A) =VX ( 12 ) 

CUR (3 , I A) =VX ( 13 ) 

CUR (4 , IA) =VX ( 14 ) 

1 CONTINUE 

c (Carry out actual integral . . . . ) 

DO 3 IFF=1 , 4 

IF (CUR(IFF,1) .EQ.O.O) OFLX ( IFF) =-999 
IF (CUR (IFF, 1) .EQ.O.O) GO TO 3 
DO 4 IA=1 , N 

4 CUR1 (IA) =CUR (IFF, IA) 

CALL INTEG ( CUR1 , XX , N , OF ) 

OFLX (IFF) =ALOG10 (OF*4 . *3 . 14159) 

3 CONTINUE 
RETURN 

c As stated this is the section when L is outside range 

55 DO 56 1=1,4 

56 OFLX ( I ) =-999 
RETURN 
END 

SUBROUTINE SRADBL (B, Bcc, Bee , AL, ALAT, RS , WLONG, ALPHA, E) 

PROGRAM TO FIND HIGH ENERGY FLUXES OF ELECTRONS AND PROTONS 
FIE=I (CM**2-S-SR) **-l FOR ELECTRONS 
FDE=DI/DE (CM**2-S-SR-MEV) **-l FOR ELECTRONS 
FIP=I (CM**2-S-SR) **-l FOR PROTONS 
FDP=DI/DE (CM**2-S-SR-MEV) **-l FOR PROTONS 

COMMON /RADAT/Brad, Btheta , Bphi , Bt , AL1 , ALATC , 

1 BE , BC , BM, THETA, FIE , FDE , FIP , FDP 

C** Corrected 09/2000 (more accurate) HBG 
RDD=. 017453292 

be=bee 

bc=bcc 

c Is ALPHA along the field line? If so, flux = 0. 

IF (ALPHA. EQ.O.O) GO TO 49 
IF ( ALPHA. EQ. 180. ) GO TO 49 

C Is BM<=BC? If not, then flux = 0. 

BM=B/ (SIN (ALPHA*RDD) **2) 

IF (BM.GT.BC) GO TO 49 

C FIND ELECTRON INTEGRAL /DIFFERENTIAL INTENSITIES 
call elntensity (alpha, E, Be, B, AL, FIE, FDE) 
call plntensity (alpha, E , Be, B, AL, FIP, FDP) 

GO TO 50 

c Zero flux trap. . . . 

49 FIE=0 
FDE=0 
FIP=0 
FDP=0 
BM=999 . 

50 CONTINUE 

RETURN 

END 


Subroutine elntensity (pitch, Energy, Beqmag , Bmag, Lvalue , 
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c 

c 


% eiintst , edintst ) 

Calculate the Intensity for a specified pitch angle 
and (r , col at , elong) geographic position, 
real E2x, E3x, Energy, pitch, AA(4) , Lvalue 

c Get A coefficients: 

eiintst=0 . 0 
edintst=0 . 0 

call eAcoef f (pitch, AA, E2x, E3x, Beqmag, Bmag, Lvalue, iflagout ) 
if (iflagout .ne . 0) return 

xl= AA ( 1 ) - AA ( 2 ) *alogl0 (Energy) 

x2= . 5* (AA(2) - AA ( 3 ) ) *alogl0 ( (Energy* *2+E2x**2 ) / (1 . + E2x**2) ) 
x3= (AA ( 3 ) -AA ( 4 ) ) *alogl0 ( ( Energy* *3+E3x* *3 ) / (1 .+E3x**3) ) /3 . 
eiintst = 10 . ** (xl+x2+x3) 

edintst= (eiintst /Energy) * (AA(2) + (AA (3 ) -AA (2) ) / 

% (1 .+ (E2x/Energy) **2) + (AA(4) -AA(3) )/(!.+ (E3x/Energy) **3) ) 


Return 

End 


Subroutine eAcoef f (pitch 
real Lvalue 

real AA(4), An(4,27), Bn 
real E2 (27) , E3 (27) , E2x 
data (An (i , 1) , 
data (An (i , 2) , 
data (An (i , 3) , 
data (An (i, 4) , 

(An (i , 5) , 

(An (i , 6) , 

(An (i , 7) , 
data (An (i , 8) , 
data (An (i , 9) , 
data (An (i , 10) , 

(An (i , 11 ) , 

(An ( i , 12 ) , 
data (An (i , 13 ) , 
data (An (i , 14 ) , 

(An(i, 15) , 

(An (i , 16 ) , 

(An ( i , 17 ) , 
data (An (i , 18 ) , 
data (An (i , 19) , 
data (An (i , 20) , 
data (An (i , 21 ) , 
data (An (i , 22 ) , 

(An(i, 23) , 

(An ( i , 24 ) , 

(An(i, 25) , 

(An (i , 26 ) , 

(An(i, 27) , 


, AA, E2x, E3x, Beqmag, Bmag, Lvalue, iflagout) 

(4,27), pitch, Ax (4 ) , Bx(4) 

, E3x, Ltable (27) 


data 

data 

data 


data 

data 


data 

data 

data 


data 

data 

data 

data 

data 


data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 

data 


(Bn (i , 1) , 
(Bn (i , 2 ) , 
(Bn (i , 3) , 
(Bn (i , 4 ) , 
(Bn (i , 5) , 
(Bn (i , 6) , 
(Bn (i , 7) , 
(Bn (i , 8) , 
(Bn (i , 9) , 
(Bn(i, 10) 
(Bn(i, 11) 
(Bn(i, 12) 
(Bn(i, 13) 
(Bn (i , 14 ) 
(Bn(i, 15) 
(Bn (i , 16 ) 
(Bn (i , 17 ) 
(Bn (i , 18 ) 


i=l, 4) 

/ 4. 

74, 

i=l, 4) 

/5. 

02, 

i=l, 4) 

/ 4- 

91, 

i=l, 4) 

/ 4- 

98, 

i=l, 4) 

/5. 

oo, 

1-1,4) 

/ 5 . 

49, 

i=l, 4) 

/5. 

56, 

1=1,4) 

/ 5. 

51, 

i=l, 4) 

/5. 

52, 

i=l , 4 ] 

1/5. 

46, 

i=l , 4 ! 

1/5. 

38, 

1=1,4] 

1/5. 

31, 

i=l, 41 

1/5. 

.03, 

i=i , 4 ; 

1/5. 

50, 

i=l , 4 1 

1/4. 

.88, 

i-1,41 

1/4. 

.67, 

1-1.41 

1/4. 

■74, 

i-1,41 

>/3. 

.58, 

i-1,41 

1/3. 

.50, 

i-1,41 

1/4. 

.57, 

i-1,41 

>/3. 

.26, 

i-1,41 

)/3. 

.95, 

i-1,41 

>/3. 

.33, 

i-1,41 

)/3. 

.42, 

i-1,41 

>/3. 

.02, 

i-1,41 

)/3. 

.64, 

i-1,41 

)/2. 

.09, 

i=l, 4) 

/ 4 

.74, 

1=1,4) 

/ 5 

.02, 

1=1,4) 

/ 4 

■91, 

i=l, 4) 

/ 4 

.98, 

i=l, 4) 

/ 5 

-00, 

i=l, 4) 

/ 5 

•49, 

i=l, 4) 

/ 5 

. 56, 

i=l, 4) 

/ 5 

•51, 

i=l, 4) 

/ 5 

. 52, 

i=l , 4 

) / 5 

.46, 

i=l,4 

) / 5 

-38, 

i-1.4 

) / 5 

.48, 

i-1,4 

) / 5 

.74, 

i-1,4 

) / 4 

. 90, 

i=l, 4 

) / 4 

•72, 

i=l , 4 

) / 4 

•50, 

i=l , 4 

) / 4 

•12, 

i-1,4 

) / 5 

.30, 


0.00, 

0.00, 

6.47/ 

0.00, 

0.00, 

4.01/ 

0.00, 0.00, 

4.86/ 

0.00, 0.00, 

4.41/ 

0.00, 

0.00, 

5.37/ 

0.00, 

0.23, 

3.00/ 

0.00, 

0.34, 

3.00/ 

0.00, 0.27, 

4.48/ 

0.00, 0.25,. 

3 . 11/ 

0.00, 

0.28, 

3.42/ 

0.00, 

0.27, 

4 . 96/ 

0.00, 

0.29, 

5.23/ 

0.00, 0.20, 

5.48/ 

0.00, 0.00, 

4.27/ 

0.00, 

0.20, 

5.71/ 

0.00, 

0.00, 

4.29/ 

0.00, 

0.00, 

4.23/ 

0.00, 0.34, 

5.74/ 

0.00, 

0.27, 

5.26/ 

0.00, 

0.05, 

6.67/ 

0.00, 

0.00, 

4.22/ 

0.00, 

0.06, 

3 . 00/ 

0.00, 

0.00, 

6.10/ 

0.00, 

0.00, 

5.50/ 

0.00, 

0.83, 

6 . 00/ 

0.00, 

1.17, 

3.01/ 

2.23, 

0.00, 

7.59/ 

0.00, 

0.00, 

6.47/ 

0.00, 

0.00, 

4 . 01/ 

0.00, 

0.00, 

4.86/ 

0.00, 

0.00, 

4.41/ 

0.00, 0.00, 

5.37/ 

0.00, 

0.23, 

3.00/ 

0.00, 

0.34, 

3 . 00/ 

0.00, 

0.27, 

4.48/ 

0.00, 

0.25, 

3.11/ 

0.00, 0.28, 

3.42/ 

0.00, 

0.27, 

4.96/ 

0.00, 

0.29, 

5.23/ 

0.00, 

0.20, 

5.48/ 

0.00, 

0.61, 

4.27/ 

0.21, 0.78, 

5.71/ 

0.00, 

0.56, 

4 . 29/ 

0.00, 

0.68, 

4.23/ 

0.00, 0.11, 

5.74/ 
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500 

501 


51 

510 


c 


520 


c 

c 


data 

(Bn (i , 19) , i=l , 4 ) /5 . 38 , 

0.00 

, 0.11, 

5.26/ 

data 

(Bn (i , 20) , i=l,4)/3.73. 

1.01 

, 0.68, 

6.67/ 

data 

(Bn(i, 21) , i=l,4)/4.60. 

0.26 

, 0.00, 

4.22/ 

data 

(Bn(i, 22) , 1=1,4)74.83, 

0.00 

, 0.19, 

3 . 00/ 

data 

(Bn(i,23) , i=l , 4 ) /3 . 20 , 

0.21 

, 0.66, 

6.10/ 

data 

(Bn(i, 24) , i = l , 4 ) /2 . 56 , 

0.16 

, 2.45, 

5.50/ 

data 

(Bn(i,25) , i=l , 4 ) /2 . 11 , 

0.00 

, 2.22, 

6.00/ 

data 

(Bn(i, 26) , i=l , 4 ) /I . 98 , 

0.00 

, 3-78, 

3.01/ 

data 

(Bn (i , 27) , i=l , 4 ) /2 . 09 , 

2.23 

, 0.00, 

7.59/ 

data 

Ltable/2 .30,2.33,2.34,2 

.36,2 

.38,2.46,2.56, 

> 

2.93,3.04,3.08,3.20,3. 

40,3 . 

90,4.30 

,4.44,4 

r 

6.20,6.50,7.20, 8.00,9. 

50,11 

.00,13. 

00/ 

data 

E2/.6, .6, .6, .6, .6, .6, .6 

, .6, . 

2, .2, .2 

, -2, .2, 

r 

.1, .1, .1, .1,0.08,0.08, 

0.08, 

0.08,0. 

08,0.08 

data 

E3/1. 5,1. 5, 1.5, 1.5, 1.5, 

2. ,2. 

/ 2 . , 2 . , 

2. ,2. ,2 


2.,2.,1.,1.,1.,1.,1.,0.5,0.5,0.5,0.5,0.5,0.5,0.5/ 


iflagout=0 

if (Lvalue .gt .Lt able (27) .or .Lvalue .It .Ltable (1) ) then 
iflagout = 1 
goto 520 
end if 

Interpolate coeffs. 
if (Lvalue .GE. Ltable(27)) then 
nmin=26 
goto 501 
end if 

if (Lvalue .LT. Ltable (7)) then 
nmin=l 
end if 

if ( Lvalue. GE. Ltable (7) . and. Lvalue . LT. Ltable (13) ) then 
nmin=7 
end if 

if (Lvalue. GE. Ltable (13) . and. Lvalue .LT. Ltable (20) ) then 
nmin=13 
end if 

i f ( Lvalue . GE . Ltable (20). and . Lvalue . LE . Ltable (27)) then 
nmin=20 
end if 

do 500, j=nmin,26 

if (Lvalue .LT. Ltable (j+1) ) then 
nmin=j 
Goto 501 
end if 
continue 
continue 

scale= (Lvalue -Ltable (nmin) ) / (Ltable (nmin+1) -Ltable (nmin) ) 
do 510, 1=1,4 

Ax ( I ) = An (I, nmin) + scale* (An (I , nmin+1 ) -An (I , nmin) ) 

Bx ( I ) = Bn (I, nmin) + scale* (Bn (I , nmin+1 ) -Bn (I , nmin) ) 
continue 
continue 

E2x= E2 (nmin) + scale* (E2 (nmin+1) -E2 (nmin) ) 

E3x= E3 (nmin) + scale* (E3 (nmin+1) -E3 (nmin) ) 

Calculate AA 
rp=pitch *3.14159/180. 
do 520, J=1 , 4 

AA ( J) = ( (Beqmag/Bmag)* (Ax(J) - Bx ( J) ) * (sin (rp) ) **2 . ) +Bx ( J) 
continue 
Return 
End 

Subroutine plntensity (pitch. Energy, Beqmag, Bmag, Lvalue, 

% piint st , pdintst ) 

Calculate the log of the Intensity for a specified pitch angle 
and (r, lat) geographic position, 
real E2x, E3x, Energy, pitch, AA(4), Lvalue 
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c 


Get A coefficients: 
piintst=0 . 0 
pdintst=0 . 0 

call pAcoef f (pitch, AA, E2x, E3x, Beqmag, Bmag, Lvalue, iflagout) 
i f ( i f lagout . ne . 0 ) return 
xl = AA(1) - AA (2 ) *aloglO (Energy) 

x2=0 . 5* (AA ( 2 ) -AA ( 3 ) ) *aloglO ( (Energy** 2 + E2x**2)/(1 + E2x**2) ) 
x3= (1-/3-)* (AA(3) -AA(4) ) *aloglO ( ( Energy* *3+E3x* *3 ) / (1 + E3x**3) ) 
piintst =10 . ** (xl+x2+x3 ) 

pdintst= (piintst /Energy) * (AA(2) + (AA(3) -AA(2) ) / 

% (1.+ (E2x/Energy) **2) + (AA(4) -AA(3) ) / (1 .+ (E3x/Energy) **3) ) 

Return 

End 


Subroutine pAcoef f (pitch, AA, E2x, E3x, Beqmag, Bmag, Lvalue, iflagout) 
real Lvalue 

real AA (4 ) , An(4,29), Bn(4,29), pitch, Ax(4), Bx(4) 
real E2 (29) , E3 (29) , E2x, E3x, Ltable(29) 


data 

(An (i , 1) , 

i=l , 4 ) /1.86, 

0.00, 

0.00, 

4.57/ 

data 

(An(i, 2) , 

i=l, 4) /2.07, 

0.00, 

0.00, 

4 . 53/ 

data 

(An(i, 3) , 

i=l, 4) /2.17, 

0.00, 

0.00, 

6 . 09/ 

data 

(An(i, 4) , 

i=l, 4) /3.18, 

0.05, 

0.00, 

3.05/ 

data 

(An (i , 5) , 

i=l, 4) /2.14, 

1.85, 

0.00, 

3.00/ 

data 

(An (i , 6) , 

i=l, 4) /3.80, 

0.75, 

0.00, 

3.00/ 

data 

(An(i, 7) , 

i=l, 4) /3 . 37 , 

1.28, 

0.00, 

3.00/ 

data 

(An(i, 8) , 

i=l, 4) /3 . 26 , 

2.63, 

0.00, 

3 . 00/ 

data 

(An (i , 9) , 

i=l , 4 ) /3.14, 

2.93, 

0.00, 

3.00/ 

data 

(Anti, 10) 

, i=l , 4 ) /0 . 37 , 

0.00, 

0.00, 

3.00/ 

data 

(An(i, 11) 

, i=l,4)/2.91, 

0.29, 

0.00, 

3.00/ 

data 

(An (i , 12 ) 

, i=l , 4 ) /2 . 86 , 

0.02, 

0.02, 

3.00/ 

data 

(An(i, 13) 

, i=l , 4 ) /2 . 80 , 

0.98, 

0.00, 

3.00/ 

data 

(An (i , 14) 

, i = l , 4 ) /2 . 16 , 

1.77, 

0.00, 

3.00/ 

data 

(An(i, 15) 

, i = l , 4 ) / -0 . 13 , 

r 0.00, 

, 0.26, 

3.00/ 

data 

(An (i , 16 ) 

, i=l , 4 ) /I . 29 , 

0.00, 

0.00, 

3.00/ 

data 

(An(i, 17) 

, i=l , 4 ) /0 . 63 , 

0.46, 

0.29, 

3.22/ 

data 

(An(i, 18) 

, i=l , 4 ) /0 . 61 , 

5.05, 

0.20, 

3.12/ 

data 

(An(i, 19) 

, i=l,4)/0.47. 

5.25, 

0.21, 

3.55/ 

data 

(An(i, 20) . 

, i=l,4)/0.68, 

1.45, 

0.00, 

4.41/ 

data 

(An(i, 21) , 

, i=l , 4) /I . 77, 

0.00, 

0.00, 

3.00/ 

data 

(An(i, 22) , 

, i=l , 4 ) /2 . 32 , 

3.32, 

0.00, 

3.08/ 

data 

(And, 23) , 

, i=l, 4 ) /2 . 40, 

3.34, 

0.00, 

3 . 00/ 

data 

(Anti, 24) , 

, i = l , 4 ) /2 . 46 , 

2.59, 

0.52, 

4.72/ 

data 

(And, 25) , 

, i=l , 4 ) /2 . 46 , 

3.55, 

0.25, 

3 . 85/ 

data 

(And, 26) , 

, i=l, 41/2.31, 

6.28, 

1.73, 

3.35/ 

data 

(And, 27) , 

, i=l, 4)/1.90, 

6.52, 

1.81, 

4.11/ 

data 

(Anti, 28) , 

, i = l , 4 ) /1 . 24 , 

0.00, 

6.35, 

3.94/ 

data 

(Anti, 29) , 

, i=l, 4) /I . 22 , 

0.00, 

6.72, 

3.31/ 

data 

(Bn(i, 1) , 

i = l, 4) /1 . 86 , 

0.00, 

0.00, 

4.57/ 

data 

(Bnd, 2) , 

i=l , 4 ) /2.07, 

0.00, 

0.00, 

4.53/ 

data 

(Bnd, 3) , 

i=l, 4) /2.17, 

0.00, 

0.00, 

6 . 09/ 

data 

(Bn (i , 4) , 

i=l, 4) /3.18, 

0.05, 

0.00, 

3.05/ 

data 

(Bn (i , 5) , 

i=l, 4) /2 . 14 , 

1.85, 

0.00, 

3.00/ 

data 

(Bn(i, 6) , 

i=l, 4) / 3.80, 

0.75, 

0.00, 

3.00/ 

data 

(Bn (i , 7) , 

i=l, 4) /3 . 37 , 

1.28, 

0.00, 

3.00/ 

data 

(Bnd, 8) , 

i = l, 4) /3 . 26 , 

2.63, 

0.00, 

3 . 00/ 

data 

( Bn ( i , 9) , 

i=l, 4) /3.14, 

2.93, 

0.00, 

3 . 00/ 

data 

(Bn(i, 10) , 

, i=l , 4 ) /0 . 67 , 

4.87, 

6.27, 

3.00/ 

data 

(Bn (i , 11 ) , 

, i=l,4)/-0.25, 

, 6.63, 

0.00, 

3.00/ 

data 

(Bn(i, 12) , 

, i=l , 4 ) /0 . 29 , 

8.30, 

0.00, 

3 . 00/ 

data 

(Bnd, 13) , 

, i=l, 4)/0.46, 

7.41, 

0.00, 

3.00/ 

data 

(Bn (i , 14) , 

, i=l , 4 ) /0 . 40 , 

7.10, 

0.00, 

3.00/ 

data 

(Bnd, 15) , 

, i=l , 4 ) /0 . 58 , 

7.70, 

0.52, 

3.00/ 

data 

(Bn(i, 16) , 

, i=l , 4 ) /-0 . 10, 

, 9.00, 

0.15, 

3.00/ 

data 

(Bn(i, 17) , 

, i=l , 4 ) /0 . 16 , 

8.55, 

o.oo, 

3.22/ 

data 

(Bn(i, 18) , 

, i=l , 4 ) / -0 . 03 , 

- 6.93, 

0.00, 

3.12/ 

data 

(Bn(i, 19) , 

, i=l , 4 ) / -0 . 06 , 

- 7.15, 

0.00, 

3.55/ 

data 

(Bnd, 20) , 

, i=l , 4 ) /-I • 13 , 

- 6.73, 

6.85, 

4 .41/ 

data 

(Bnd, 21) , 

, i=l, 4)7-0.98, 

i 6.78, 

8.25, 

3.00/ 
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data 

(Bn(i,22) , 

i=l,4)/0.10. 

5.51, 

2.99, 

3.08/ 

data 

(Bn(i,23) , 

i=l,4)/0.06. 

5.81, 

2.66, 3.00/ 

data 

(Bn (i , 24 ) , 

i=l,4)/0.30, 

5.83, 

0.12, 4.72/ 

data 

(Bn(i,25) , 

i = l < 4 ) /0 . 81 , 

4.90, 

0.61, 

3.85/ 

data 

(Bn (i , 26 ) , 

i=l,4) /1.07, 

0.48, 

1.73, 

3.35/ 

data 

(Bn (i , 27) , 

i=l,4) /0.74, 

0.00, 

0.85, 4.11/ 

data 

(Bn(i, 28) , 

i=l , 4 ) /I . 24 , 

0.00, 

6.35, 

3.94/ 

data 

(Bn (i , 29) , 

i=l , 4 ) /I . 22 , 

0.00, 

6.72, 

3.31/ 


data Ltable/2. 30 ,2.32 ,2. 34 ,2. 43 ,2. 52 ,2. 67 ,2. 77 ,2. 82 ,2. 86, 

* 3.10,3.20,3.31,3.44,3.71,3.99,4.13,4.27,4.38,4.49,5.01,5.36, 

* 6.20,6.45,6.70,7.63,8.40,10.30,12.00,13.00/ 

data E2/1. , 1. ,1.,1.,1.,1.,0.5,0.3,0.3,1.,1.,0.5,0.5,0.5,0.5,0.5, 

* 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.3, 0.3, 0.3, 0.3/ 
data E3/60., 60. ,60. ,60. ,60. ,60. ,60. ,60. ,60. ,60. ,60. ,60. ,60., 

* 60. ,60. ,60. ,60. ,60. ,60. ,0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8/ 

iflagout=0 

if (Lvalue. gt.Ltable (29) .or. Lvalue .It. Ltable(l)) then 
iflagout = 1 
goto 520 
end if 

c Interpolate coeffs. 

if (Lvalue .GE. Ltable(29)) then 
nmin=28 
goto 501 
end if 

if (Lvalue .LT. Ltable(7)) then 
nmin=l 
end i f 

if ( Lvalue. GE.Lt able (7) . and. Lvalue . LT. Lt able (13) ) then 
nmin=7 
end if 

if (Lvalue .GE.Lt able (13) . and. Lvalue. LT.Ltable (20) ) then 
nmin=13 
end if 

if (Lvalue .GE . Ltable (20) . and. Lvalue .LE .Ltable (29) ) then 
nmin=20 
end if 

do 500, j=nmin,28 

if (Lvalue .LT. Ltable (j+1) ) then 
nmin= j 
Goto 501 
end if 

500 continue 

501 continue 

scale= (Lvalue -Ltable (nmin) ) / (Ltable (nmin+1) -Ltable (nmin) ) 
do 510, 1=1,4 

Ax ( I ) = An (I, nmin) + scale* (An (I , nmin+1) -An (I , nmin) ) 

Bx ( I ) = Bn (I, nmin) + scale* (Bn (I, nmin+1) -Bn (I, nmin) ) 

51 continue 

510 continue 

E2x= E2 (nmin) + scale* (E2 (nmin+1) -E2 (nmin) ) 

E3x= E3 (nmin) + scale* (E3 (nmin+1) -E3 (nmin) ) 
c Calculate AA 

rp=pitch *3.14159/180 
do 520, J=1 , 4 

AA(J) = ( (Beqmag/Bmag) * (Ax( J) - Bx ( J) ) * (sin (rp) ) **2 ) +Bx ( J) 

520 continue 

Return 
End 

SUBROUTINE INTEG ( YY , XX, N, AINT) 

DIMENSION YY ( 1 ) , XX ( 1 ) 

AINT =0 
DO 1 1=1 , N 
1 AINT= YY ( I ) +AINT 

AINT= (AINT- ( YY ( 1 ) +YY (N) ) * .5) * (XX (2) -XX (1) ) 

RETURN 

END 
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Input file — “traj.in” 

The input file provides the time, trajectory (rs = radial distance in R s , alat = 
latitude in degrees, wlong = west longitude in degrees), and magnetic field components 
(ali = L-shell, bsc = magnetic field amplitude in gauss at spacecraft, beq = magnetic field 
amplitude in gauss at magnetic equator, befit = critical or cut-off magnetic field 
amplitude at top of atmosphere). “51” is the number of points and “60” is the time step 
in seconds required to calculate the fluence. The last two rows are for program control. 

Test Case 2 for Saturn — equator 


51 60 


time 

rs 

alat 

wlong 

ali 

bsc 

beq 

berit 

bm 

60.00000 

2.30000 

.00000 

.00000 

2.30000 

.01708 

.01706 

.43054 

.21536 

120.00000 

2.32000 

.00000 

.00000 

2.32000 

.01666 

.01663 

.43272 

.21536 

180.00000 

2.33000 

.00000 

.00000 

2.33000 

.01645 

.01643 

.43381 

.21536 

240.00000 

2.34000 

.00000 

.00000 

2.34000 

.01624 

.01622 

.43488 

.21536 

300.00000 

2.36000 

.00000 

.00000 

2.36000 

.01584 

.01582 

.43701 

.21536 

360.00000 

2.38000 

.00000 

.00000 

2.38000 

.01545 

.01544 

.43911 

.21536 

420.00000 

2.43000 

.00000 

.00000 

2.43000 

.01454 

.01452 

.44423 

.21536 

480.00000 

2.46000 

.00000 

.00000 

2.46000 

.01403 

.01401 

.44723 

.21536 

540.00000 

2.52000 

.00000 

.00000 

2.52000 

.01307 

.01305 

.45305 

.21536 

600.00000 

2.56000 

.00000 

.00000 

2.56000 

.01247 

.01246 

.45680 

.21536 

660.00000 

2.67000 

.00000 

.00000 

2.67000 

.01102 

.01101 

.46666 

.21536 

720.00000 

2.70000 

.00000 

.00000 

2.70000 

.01066 

.01065 

.46923 

.21536 

780.00000 

2.77000 

.00000 

.00000 

2.77000 

.00989 

.00988 

.47505 

.21536 

840.00000 

2.82000 

.00000 

.00000 

2.82000 

.00938 

.00937 

.47906 

.21536 

900.00000 

2.83000 

.00000 

.00000 

2.83000 

.00928 

.00928 

.47985 

.21536 

960.00000 

2.86000 

.00000 

.00000 

2.86000 

.00900 

.00899 

.48218 

.21536 

1020.00000 

2.93000 

.00000 

.00000 

2.93000 

.00838 

.00837 

.48747 

.21536 

1080.00000 

3.04000 

.00000 

.00000 

3.04000 

.00751 

.00751 

.49536 

.21536 

1140.00000 

3.08000 

.00000 

.00000 

3.08000 

.00723 

.00722 

.49811 

.21536 

1200.00000 

3.10000 

.00000 

.00000 

3.10000 

.00709 

.00709 

.49946 

.21536 

1260.00000 

3.20000 

.00000 

.00000 

3.20000 

.00645 

.00645 

.50599 

.21536 

1320.00000 

3.31000 

.00000 

.00000 

3.31000 

.00584 

.00583 

.51277 

.21536 

1380.00000 

3.40000 

.00000 

.00000 

3.40000 

.00539 

.00539 

.51802 

.21536 

1440.00000 

3.44000 

.00000 

.00000 

3.44000 

.00521 

.00520 

.52028 

.21536 

1500.00000 

3.71000 

.00000 

.00000 

3.71000 

.00416 

.00416 

.53432 

.21536 

1560.00000 

3.90000 

.00000 

.00000 

3.90000 

.00359 

.00358 

.54310 

.21536 

1620.00000 

3.99000 

.00000 

.00000 

3.99000 

.00335 

.00335 

.54697 

.21536 

1680.00000 

4.13000 

.00000 

.00000 

4.13000 

.00302 

.00302 

.55268 

.21536 

1740.00000 

4.27000 

.00000 

.00000 

4.27000 

.00274 

.00274 

.55800 

.21536 

1800.00000 

4.30000 

.00000 

.00000 

4.30000 

.00268 

.00268 

.55910 

.21536 

1860.00000 

4.38000 

.00000 

.00000 

4.38000 

.00254 

.00254 

.56195 

.21536 

1920.00000 

4.44000 

.00000 

.00000 

4.44000 

.00244 

.00244 

.56402 

.•21536 

1980.00000 

4.49000 

.00000 

.00000 

4.49000 

.00236 

.00236 

.56569 

.21536 

2040.00000 

4.96000 

.00000 

.00000 

4.96000 

.00175 

.00175 

.57971 

.21536 

2100.00000 

5.01000 

.00000 

.00000 

5.01000 

.00170 

.00170 

.58103 

.21536 

2160.00000 

5.36000 

.00000 

.00000 

5.36000 

.00139 

.00139 

.58951 

.21536 

2220.00000 

5.41000 

.00000 

.00000 

5.41000 

.00135 

.00135 

.59061 

.21536 

2280.00000 

5.87000 

.00000 

.00000 

5.87000 

.00106 

.00106 

.59972 

.21536 

2340.00000 

6.20000 

.00000 

.00000 

6.20000 

.00090 

.00090 

.60521 

.21536 

2400.00000 

6.45000 

.00000 

.00000 

6.45000 

.00080 

.00080 

.60887 

.21536 

2460.00000 

6.50000 

.00000 

.00000 

6.50000 

.00078 

.00078 

.60956 

.21536 

2520.00000 

6.70000 

.00000 

.00000 

6.70000 

.00071 

.00071 

.61215 

.21536 

2580.00000 

7.20000 

.00000 

.00000 

7.20000 

.00057 

.00057 

.61768 

.21536 

2640.00000 

7.63000 

.00000 

.00000 

7.63000 

.00048 

.00048 

.62152 

.21536 

2700.00000 

8.00000 

.00000 

.00000 

8.00000 

.00042 

.00042 

.62423 

.21536 

2760.00000 

8.40000 

.00000 

.00000 

8.40000 

.00036 

.00036 

. 62664 

.21536 

2820.00000 

9.50000 

.00000 

.00000 

9.50000 

.00025 

.00025 

.63108 

.21536 

2880.00000 

10.30000 

.00000 

.00000 

10.30000 

.00020 

.00020 

.63598 

.21536 

2940.00000 

11.00000 

.00000 

.00000 

11.00000 

.00016 

.00016 

.63975 

.21536 

3000.00000 

12.00000 

.00000 

.00000 

12.00000 

.00012 

.00012 

.64447 

.21536 

3060.00000 

13.00000 

.00000 

.00000 

13.00000 

.00010 

.00010 

.64852 

.21536 

1 

0 

0 

0 

0 

0 

0 

0 


2 

0 

0 

0 

0 

0 

0 

0 
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Sample Output Data — Sfluence.dat 

Electron and proton integral and differential fluences at specified energies. 

Divine IOM 5217-90-029, 02/28/90 
FLUENCE SPECTRA 

Electrons Protons 


Energy 

Integral 

Differential 

Integral 

Differential 

(MeV) 

(cm**-2 ) 

(cm* *2 MeV) **-l 

(cm**-2) 

(cm**2 MeV)** 

. 10 

0 . 987E+10 

0 . 143E+11 

0 . 315E+10 

0 . 157E+12 

.20 

0 . 882E+10 

0. 837E+10 

0 . 191E+09 

0 . 291E+10 

.30 

0 . 806E+10 

0 . 702E+10 

0 . 685E+08 

0 . 4 68E+09 

.50 

0 . 682E+10 

0 . 541E+10 

0 . 309E+08 

0 . 701E+08 

1.00 

0 . 463E+10 

0 . 380E+10 

0 . 183E+08 

0 .807E+07 

2.00 

0 . 175E+10 

0 . 188E+10 

0 . 154E+08 

0 . 945E+06 

3.00 

0 . 616E+09 

0 . 603E+09 

0.149E+08 

0 .278E+06 

5.00 

0 . 123E+09 

0 . 804E+08 

0 . 14 6E+08 

0 . 665E+05 

10.00 

0 . 126E+08 

0 . 410E+07 

0 . 144E+08 

0 . 295E+05 

30.00 

0 . 382E+06 

0 . 400E+05 

0 . 128E+08 

0 . 148E+06 

100.00 

0 . 908E+04 

0 . 280E+03 

0.245E+07 

0 . 613E+05 

300.00 

0 . 310E+03 

0 . 317E + 01 

0 . 107E+06 

0 .107E+04 

1000.00 

0 . 783E+01 

0 . 239E-01 

0 . 289E+04 

0 . 869E+01 
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Sample Output Data — SOrbitfluence.dat 

Electron and proton integral fluences for input into a typical radiation program 

Divine, IOM 5217-90-029, 02/28/90 


FLUENCE SPECTRA 
3060.00 

Electrons 


Energy 

Integral 

(MeV) 

(cm** -2) 

.10 

0 . 987E+10 

.20 

0 . 882E+10 

.30 

0 . 806E+10 

.50 

0 . 682E+10 

1.00 

0 .463E + 10 

2.00 

0 . 175E+10 

3.00 

0 . 616E+09 

5.00 

0 . 123E+09 

10.00 

0 . 126E+08 

30.00 

0 . 382E+06 

100.00 

0 . 908E+04 

300.00 

0 . 310E+03 

1000.00 

0 . 783E+01 

Protons 

Energy 

Integral 

(MeV) 

(cm** -2) 

.10 

0 . 315E + 10 

.20 

0 . 191E+09 

.30 

0 . 685E+08 

.50 

0 . 309E+08 

1.00 

0 . 183E+08 

2.00 

0 . 154E+08 

3.00 

0 . 149E+08 

5.00 

0 . 146E+08 

10.00 

0 . 144E+08 

30.00 

0 . 128E+08 

100.00 

0 .245E + 07 

300.00 

0 . 107E+06 

1000.00 

0.289E+04 
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Sample Output Data — SfluxdifT.dat 
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iPPENDIX VI. Original Divine Model Report 


JET PROPULSION LABORATORY INTEROFFICE MEMORANDUM 

5217-90-029 

1990 February 28 

TO: A. R. Hoffman (301 -466) 

FROM: Neil Divine (301 -460) 

SUBJECT: Numerical Models for Electron and Proton Distributions in Saturn’s Radiation 

Belts 

SUMMARY 


Published data from charged particle experiments aboard three spacecraft during their 
flybys at Saturn have been used to generate numerical models for the radiation belts 
between 2.3 and 13 Saturn radii. The models describe the electron distributions at 
energies between 0.04 and 10 MeV and the proton distributions at energies between 0.14 
and 80 MeV. They are intended to predict particle intensity, flux and fluence for the 
Cassini orbiter. Review comments on the models will be solicited. 
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1. INTRODUCTION 

An extensive data set describing the distributions of energetic electrons and ions in 
Saturn’s radiation belts has become available following three spacecraft flybys in 1979, 
1980, and 1981. These data are presented in numerous articles in the research literature 
(see Sec. 2 below) and are summarized by Van Allen, 1984 (reference VA84, introducing the 
citation scheme employed in this memorandum). The following sections describe numerical 
models for the distributions based on a selection of published profiles within the data set. 
The present models are intended to replace those derived from rather limited, and to some 
extent preliminary, subsets of the data, as described in DI80, EV87a, and EV87b. 

Table 1 specifies values for several Saturn quantities used with the models. The five 
magnetic field coefficients represent the zonal harmonic model (Z 3 ); all other coefficients are 
set equal to zero in this description of the magnetic field from C084 (based on numerical fits 
to magnetic field data from the flyby spacecraft). In the equations (1 through 3) of this 
reference, a is equal to 1.0 R s (Table 1), and the longitudes <t> and y are both measured 
positive east from the prime meridian of the planet, including prograde rotation with rate <» 
(Table 1). Although the internal terms (proportional to g{\ g$, and g$) are axisymmetric in 
this representation, the external ones (proportional to Gp and Gj ) are not (because G/*0). 
For a given position with respect to Saturn, the magnetic shell parameter L is obtained by 
evaluating the magnetic field vector 7? at the point, tracing the field line to the magnetic 
equator (where B e is the minimum value of |Z?| along the line) and beyond to a location 
where |i?| has the same value as at the initial point, evaluating the adiabatic invariant 7 as an 
integral along the field line (RO70, p. 48, eq. 2.38), and using a series expansion 
approximation to obtain L (RO70, pp. 154-155, where k a is the same as g? in Table 1 here). 
Appendix A provides additional detail for the properties of the magnetic field description and 
of charged particle motion in the radiation belts. 


2. FLYBY SPACECRAFT AND THE DATA PROFILES 

Table 2 specifies Keplerian elements and other data for the trajectories of the three 
Saturn flyby spacecraft. The references cited in the last row of Table 2 are the sources for 
all the entries except the longitudes, which were estimated from Figures 2 and 3 of C084 
(pp. 358-359; note that these are west longitudes). Calculation of the flyby trajectories using 
Keplerian orbits about Saturn (as a point mass), including both planetary rotation and the 
magnetic field description above (Sec. 1 , last para.) yield satisfactory agreement (to within 
about 0.05 in L) with the magnetic coordinates (and times) shown in Figure 1, for L<9. 
Attempts to simplify this process (e.g., by omitting the external field harmonics) yielded 
significantly degraded agreement (e.g., discrepancies near 0.2 in L). The trajectory 
characterizations in the following paragraphs use the entries in Tables 1 and 2 (not those 
from the other references or from Fig. 1 ). 
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The Pioneer 11 trajectory at Saturn was nearly equatorial, and included a minimum 
L=1 .43 (maximim B=0.0801 G). The charged particle instrumentation aboard Pioneer 11 
consists of four subsystems, namely the Geiger Tube Telescope (GTT, provided by the Univ. 
of Iowa), Charged Particle Telescope (CPT, Univ. of Chicago), Trapped Radiation Detector 
(TRD, Univ. of California at San Diego), and Cosmic Ray Telescope (CRT, Goddard Space 
Flight Center and Univ. of New Hampshire). These four instruments were completed prior to 
the 1 973 launch and occupy fixed orientations aboard a spin-stabilized spacecraft intended 
primarily for data acquisition at its Jupiter flyby. 

The Voyager 1 trajectory at Saturn was highly inclined, and included about two hours at 
nearly constant L (5.9±0.1) inbound, maximum 5=0.011 G, and minimum L=4.36. The 
charged particle instrumentation aboard both Voyager 1 and 2 consists of two subsystems, 
namely the Low-Energy Charged Particle Experiment (LECP, provided by Applied Physics 
Laboratory, Johns Hopkins Univ.) and the Cosmic Ray Subsystem (CRS, Goddard Space 
Flight Center). These instruments were completed prior to the 1977 launch and occupy 
either fixed (CRS) or commandable rotating orientations (LECP) aboard three-axis-stabilized 
spacecraft intended for data acquisition at both Jupiter and Saturn. 

The Voyager 2 trajectory at Saturn was moderately inclined, and included maximum 
5=0.012 G and minimum L= 2.77. There were also pairs of magnetic coordinates (5,L) 
matched twice along the Voyager trajectories, namely (0.007, 4.43), (0.0036, 4.47), (0.0029, 
6.06), (0.0017, 6.55) and (0.001, 8.45); these points (see crossings in Fig. 1) enable checks 
for the reproducibility of the data (as in KR83, p. 8888, Fig. 17 and 18). 

Tables 3 and 4 list those profiles which have been selected from the published data 
(acquired during these three flybys) and used in the present electron and proton models 
respectively. The extensive notes accompanying these tables suggest some of the diversity 
and subtlety of the results. Appendix B includes a discussion of species identification for 
some of the profiles in Tables 3 and 4, whereas Appendix C includes summaries of the 
selection criteria for these profiles and of the procedures used for their processing in the 
models. 

The histograms in Figures 2 and 3 display the coverage of the selected profiles in 
magnetic field and energy threshold at each of several values for L , namely those chosen for 
use as fitting points in the models (see App. C for a discussion of the procedures used). 
Even within the ranges spanned by the profiles the distribution of the data is quite uneven, 
as restricted by the spacecraft trajectories and by the detector thresholds, dynamic ranges, 
and backgrounds. 
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3. ALGEBRAIC FORMULATION OF THE MODELS 


The basic independent variables used in these models are magnetic shell parameter L , 
local magnetic field strength B , pitch angle a, and (electron or proton) kinetic energy E . The 
form chosen to match the data (Sec. 2 and Tables 3 and 4) resembles that used in DG83 
(eq. 7 through 9) and takes the integral intensity / (in cm _2 s _1 sr _1 , for particles having 
energy >E) as its basic dependent variable. The dependence of / on E is given by 




(i) 


Differentiation of equation (1) with respect to energy leads to the differential intensity i 
(sample units cm^s^sr^MeV -1 ) in the form 


Jr ^ If . A,- A, A»-A». "| 

JE E L A » + !+(£*/£)* i + (F 3 /e; 3 4 * J 


(2) 


The dependence of / and i on a and B is given by 


Ay, = ^'(Ai^ a -) Z (a n -i) r ,) + (3) 

for subscript «= 0 , 1 , 2, and 3 (see eq. 1 and 2). The dependence of I and i on L is 
achieved by taking each of the ten parameters a n , b n , E 2 and £ 3 as linear in L between 
values specified at several fit points in L. 


Examination of equations (1) through (3) reveals the following interpretations (illustrating 
the roles of the ten model parameters) at a given value of L : 

(1 ) the spectra resemble power law segments in which / is proportional to E* 7 , where the 
integral exponent y=>-d (log/ )td (log£ ) has the value A { for kinetic energy E«E 2 , the 
value A 2 for £ 2 «E«£ 3 , and the value A 3 for E»E 3 (here y equals the quantity in 
brackets in eq. 2); 


(2) the coefficients A n vary from a n for pitch angle a=90 deg at the magnetic equator 
{B=B t ) to b H for pitch angle a=0 at any B >B e , and to b n for any pitch angle ct<90 deg at 
high magnetic latitude (S e /B=0); 


(3) at energy £=1 MeV, and at the magnetic equator (B=B e ), two parameters can be 
identified as a 0 =log I x and as 6 0 =log / (J , where subscripts distinguish intensities 
perpendicular to (ce=90 deg) and parallel to (a=0) the magnetic field and 

(4) at high magnetic latitude (B e /B =0), the intensity I is isotropic (independent of pitch angle 

a) and has the same v- je as 1^ at the equator. 
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In this formulation the kinetic energy £ and the transition energy parameters £ 2 and £ 3 all 
have units MeV, and A 0 , a g and b a are evaluated so as to yield / in the units cm _2 s -1 sr _l . 

As a subsidiary dependent variable consider the integral, omnidirectional flux J , as given 


where the last equality uses the pitch angle cosine, jx=(cosa). The corresponding 
differential, omnidirectional flux is given by 


Some additional quantities of interest which may be computed from the model are introduced 
in Appendix A. 

To understand the response of a detector in the presence of the charged particle 
distributions, consider the general expression 


Here the sum proceeds over the particle species (each having local differential intensity i , as 
in eq. 2), the solid angle integral proceeds over all directions as seen from the detector 
(subscript D ), and Q D is some quantity induced in the detector per incident particle. The 
cross-section o D may depend not only on direction of particle arrival and on its kinetic energy 
(i depends on these two as well, of course), but also on detector discrimination level. 
Analysis of the TRD data from Pioneer 1 1 (as presented in FM80) includes examples in 
which a single detector is used with three discrimination levels (channels), as well as an 
example in which Q D is charge (leading to £ as a current, as contrasted with the more 
common Q D = 1 leading to £ as a countrate). For the present model we assume that 
additional complications (not reflected in eq. 6) such as dead-time correction and 
coincidence or anti-coincidence discrimination (pp. 8933-8934 of SM83 provide a thorough 
sample discussion) are absent or already accounted for in the published data profiles (e.g., 
those in Tables 3 and 4) for £ . 


by 




( 5 ) 



( 6 ) 
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The present models employ considerable further simplification of equation (6) based on 
the approximations (in common use among the spacecraft investigators) that 

(1 ) only a single species contributes significantly to R at a given time, 

(2) average values (independent of direction) may be used for Q D and o D , regardless of 
any anisotropies in i , and 

(3) constant values may be used for Q D and o D above some discrimination threshold and 
between particle kinetic energies E a and E b , with o D =0 elsewhere. 

With these approximations equation (6) becomes 


where subscripts a and b imply evaluation of equation (4) for omni flux at the two energies 
E a and E b . Only on this basis can the conversion value RIJ in column 6 of Tables 3 and 4 
be taken as a constant, independent of the local particle distributions, for a given profile 
(detector and channel). The conversion value is dimensionless when the response R (eq. 7) 
is a flux, and has units cm 2 , sr~\ sr -1 MeV -1 , or cm 2 C whenever the response R is 
countrate, integral intensity, differential intensity, or current (examples for each may be found 
in the notes to Tables 3 and 4). 


4. NUMERICAL EVALUATION OF THE MODELS 

Table 5 specifies the parameter values adopted at 27 values of L for the electron model. 
Figures 4 and 5 illustrate the dependence on L of the electron integral intensity and omni 
flux as evaluated using the parameter values (Table 5) in equations (1), (3) and (4). Figures 
6 and 7 illustrate the dependences on pitch angle, on magnetic field strength, and on energy 
of the electron intensity and flux. 

Table 6 specifies the parameter values adopted at 29 values of L for the proton model. 
Figures 8 through 11 illustrate the proton model dependences (cf. Fig. 4 through 7 for the 
electrons). 
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The parameter values which appear in Tables 5 and 6 were obtained by an iteration 
process whose goal was to best match the model-predicted responses (defined below) to the 
data profiles (summarized in Tables 3 and 4). The iteration uses a computer program to find 
the values of several variables x, for which some single variable Y (a function of all the x,) 
has a local minimum. In an iteration, as carried out independently at each of several values 
of L (namely those in column 1 of Tables 5 and 6), the variables x, are the parameters a n , 
b n , e 2 , and £ 3 (cf. eq. 1 and 3), and the quantity Y is evaluated as described in the following 
paragraph. 


For a specific value of L, and for one particle species (i.e., electrons or protons), several 
values of the observed countrate R are obtained from measurement of the profiles in Tables 
3 or 4 (as described in App. C, para. 4). These countrates are here identified with subscript 
D , and each such value R D is associated with a pair of energies (E a and E b from column 5 
of Table 3 or 4), with a conversion constant (R/J from column 6 of Table 3 or 4), with a 
background rate (R B from column 7 of Table 3 or 4), and with a magnetic field ratio {BJB ) for 
the spacecraft trajectory and time of measurement. Corresponding to each value of R D , we 
use the model parameters (a n , b n , etc., here playing the roles of the variables x ; in the 
previous para.) and the remaining quantities (previous sentence) in equations (1), (3), (4), 
and (7) to evaluate the model-predicted rate R M (to evaluate Rm=R from eq. 7, we substitute 
the conversion RU from column 6 of Table 3 or 4 for the product Qd<Jd)- For each trio of 
rates (datum R D , model-prediction R M , and background R B ) we compute both a residual A 
and the quantity Z using 


w[o, jLjj.OWO] 

w 

vv 

► 

i 

tor 


(8) 

r a 2 

< 


IAU6. ] 

(9) 

.K 

1 

1 1 

ft) IT 




We compute the quantity Y (previous para.) as a modified RMS (root-mean-square) residual 


using 


y=[(zwz)/(zw)] ;t 


(10) 


In equation (1 0) each sum proceeds over the rate trios available at the specific value of L , 
and w is a weighting factor used for each (with very few exceptions the level weight W=1.0 
was used throughout). 
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The procedures of the foregoing two paragraphs represent the heart of the model 
development. Additional details appear in Appendix C, but we devote the remainder of this 
paragraph to a numerical example which illustrates the evaluation of the residual A (eq. 8). 
Consider a point along the Pioneer 11 trajectory inbound at L= 2.83, for which the magnetic 
field model (Table 1 and eq. 11 through 13) and the spacecraft trajectory (Table 2) yield a 
Saturnicentric distance 2.796 R s , latitude -0.27 deg, B =0.009536 G, and B e =0.009520 G at 
time -1.951 hr from perisaturn. The data were received at Earth (using Table 2) at (16:31- 
1.951 hr +86.34 min)=16:00, at which time measurement of Figure 1 of VA80c (p. 5710) 
yields a (spin-averaged) data countrate R D = 8930 s -1 for P11 GTT detector B. For this 
profile Table 3 specifies electron energies £ a =0.56 MeV and £*=21 MeV, conversion 
(£/7)=0.00127 cm 2 , and background countrate R B =0A s -1 (this estimate is from Fig. 5 of 
VA80a, p. 419, and is well below R D in our example). For the electron model, using 
parameter values from Table 5 at L= 2.83, we evaluate the four terms in equation (1) to 
obtain log /=(5.52+0.0+0.059+0.040)=5.619 at £ a =0.56 MeV and log / =(5.52+0.0-0.328- 
2.872)=2.320 at £*=21 MeV. These intensities agree with the spectrum in Figure 6 (upper 
right) and, for this isotropic distribution {b H =a n , see eq. 3 and 4), the corresponding omni 
fluxes are /=4ti/=5.22x10 6 cm~ 2 s -1 at E a =0.56 MeV and 2.62x10 s cm~ 2 s _1 at £*=21 MeV, 
which agree with the spectrum in Figure 7 (upper right). The model-predicted rate (eq. 7) is 
then R M =(R/J)(J a -J b y= 6626 s" 1 and the residual (eq. 8) is (log£ o -log£ w )=(3.95-3.82)=+0.13. 
The outbound countrate £ o =6390 s _1 (same detector, same source in VA80c) leads to a 
residual -0.02, and these two residuals are plotted as adjacent circles at log£=-0.25 and 
{B e lB)= 1.0 in Figure 12 (upper left). 

In equation (8) the first line insures that A makes no contribution to Y when both R D and 
R m are smaller than the background R B . Equation (9) insures that a peak in dZtdx, occurs 
for |A| near A a (where x, is any one of the up to ten parameters being iterated), so that the 
iteration is most effective at reducing residuals of size comparable to A a (and the hopelessly 
large residuals are deemphasized). Parameter values were accepted for inclusion in the 
model when a satisfactory minimum in Y (eq. 10) was obtained, using A„=0.3 (note that 
|A|=0.3 corresponds to a factor 2* 1 in £). 
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Many of the rows in Tables 5 and 6 demonstrate that, for each l, no more parameter 
values were iterated than could be justified by the data ranges and by their contributions to 
the reduction of the modified RMS residual Y (eq. 10). Specifically, b n =a n was enforced 
whenever it did not lead to a major increase in Y ; note that this results in the isotropy of the 
models ( b H =a n for all n=0, 1, 2, 3) for L < 3 (no data available far from the magnetic equator). 
Further, changes in E 2 and £ 3 between adjacent values of L are avoided where possible 
(partly to enable log / linear in L), and they are mostly kept within the range of the data 
threshold energies E a (see Fig. 2 and 3). Lastly a 3 and b 3 were kept >3.0, without much 
degrading the fits, to avoid unfounded extrapolation to high intensity (and flux) at large 
kinetic energy (e.g., £>10 MeV for electrons or £>100 MeV for protons) beyond the data 
coverage. 

Figures 2 and 3 suggest (at the right) the quality of the model fits, as represented by the 
residuals (eq. 8) and distinguished by magnetic shell parameter L . In each case most of the 
points from the data profiles are matched within a factor 2 ±l in R. The modified RMS 
residuals Y (eq. 10, using A c =0.3) for electrons range from 0.14 to 0.26 (corresponding 
factors for R are 1.4 to 1.8) and for protons (excepting the disaster at L= 2.52) range from 
0.11 to 0.28 (corresponding factors for/? are 1.3 to 1.9). 

Tables 7 and 8 also illustrate the quality of the model fits, this time distinguished by the 
data profiles. The range in L for which each profile is included in the model also appears in 
these tables, and the last column contains the standard RMS residual for each (retaining eq. 
8 but disregarding eq. 9 and 10). Some profiles are consistently well matched by the models 
(e.g., row 1 in Table 7 and row 10 in Table 8), whereas others are consistently badly 
matched (e.g., rows 8 and 9 in Table 8, for which the data countrate is about 10 times larger 
than the model-predicted one; see last line of eq. 8). Many of the latter "bad actors" have 
modest ranges in L , but some appear throughout the L range of the model (e.g., row 9 in 
Table 7). 

Figure 1 2 displays some sample distributions of the residuals in energy (horizontal axes), 
in magnetic field strength (vertical axes), and in value (symbol styles). The "bad actors" 
(represented by triangles) are typically well-mixed with the well-fitted points (circles), so that 
no easy fix is possible (if they were isolated the iteration could easily accommodate them, 
considering the flexibility available in eq. 1 and 3). 

The residuals cited in Tables 7 and 8 and in Figures 2, 3 and 12 suggest that a 
reasonable uncertainty estimate for the model is 2 ±l , the same as cited for the data (and/or 
its calibration) by some of the spacecraft experimenters. Some profiles consistently exceed 
this factor in comparison with the models, but the latter are based on so large a body of data 
as to be quite robust with respect to even major changes in the interpretation of a few 
profiles. 
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5. POTENTIAL IMPROVEMENTS FOR THE MODELS 

The present model is intended for review by the experiment teams responsible for the 
data cited in Table 3 and 4. Following that review it may be worthwhile to revise the model 
employing the criticisms and suggestions which are made. Thus some improvements which 
could be made now have been deferred until the later revision, as identified in the following 
paragraphs. 

Peaks in the electron model at £ =3.4, 5.87, and 6.5 (see Fig. 4) and in the proton model 
at £=2.43, 2.67, and 4.13 (see Fig. 8) do not have good resolution in £, suggesting that 
more points in £ could usefully be included in the model set (first column of Tables 5 and 6). 
Further, for electrons in the range 5<£<7, rapid changes occur between pancake (/j_ >/^ ) 
and dumbell (/j_</„ ) distributions; it would be worthwhile to check whether such is really 
demanded by the aata, and particularly whether data along the Voyager 1 trajectory (cf. 
para. 3 of Sec. 2) exerts undue influence here. 

Both models specify isotropic distributions for £< 3 (cf. para. 6 of Sec. 4), but evidence 
for strong anisotropy can be found in the spin modulation of the Pioneer 1 1 data (see for 
example FI80, Fig. 5, and VA80a, Fig. 5). Such anisotropy could probably be introduced into 
the model without increasing the residuals. In addition, for protons at £<3, a CRAND 
(cosmic ray albedo neutron decay) source suggests a strongly pancake distribution (see 
C089, p. 5-C-8, para. 3), just as occurs in the model for 3<£<4 (Fig. 8 through 11). Further 
the CRAND spectrum would have a smaller value of the power-law exponent (see y in Sec. 
3, para. 2, item 1), as suggested in FM80 (p. 5809 has y=0.8) and in VA83 (p. 6913 has 
Y=0.7 for £>120 MeV). Possibly these considerations should be included in the model (note, 
however, that y£1 leads to infinite energy density if it extends to £-»«). 

For L>10, the field line tracing (Sec. 1, para. 2, and App. A, para. 2) leads to values of B e 
which differ significantly from the ratio g ° /£ 3 . This is the source of some discrepancies in the 
B e /B coverage appearing in Figures 2 and 3 at large £ (they were calculated differently in 
different computer runs). In such cases, presumably resulting from the importance of the 
external terms in the field model, reconsideration should be given to the questions of stable 
trapping and constant intensity along the particle trajectories and drift shells. It may also be 
worthwhile to compute the relations between £ and I (adiabiatic invariant, App. A, eq. 1 7) 
using SL74 (pp. 19-25) instead of RO70 (pp. 154-155). 
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Many of the spacecraft detectors are not omnidirectional, and they respond to the local 
particle distribution with effective cross-sections (o D in Sec. 3) which vary with particle 
energy and direction (with respect to the finite angular apertures of the detectors). Thus the 
response R M which should be evaluated from the model intensity (for comparison with the 
data R d ) is more properly proportional to the integrals displayed in equation (6) than to the 
omni flux J in equation (7). Inclusion of such an improvement (to be used in the iteration for 
model parameter values) would require more detailed specifications for the detector 
calibration, for detector orientation along the three spacecraft trajectories, and/or for time 
resolution along the data profiles. To a large extent these data are indeed available 
(although not published in detail, in most cases), but their use would introduce additional 
complexity into the model development which may or may not be justifiable. 

Peak flux and fluence spectra for the Cassini orbiter have already been computed based 
on the present modeis, as reported in DI90; these results will be replaced when improved 
model distributions are completed. 
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APPENDIX A. MAGNETIC FIELD AND CHARGED PARTICLE MOTION 


The following paragraphs, which supplement the discussions in Section 1 (para. 2) and 
Section 3 (following eq. 5), represent considerable simplifications of the thorough radiation 
belt analyses presented in RO70, SL74, and C084. 


Using j? f , j? 0 , and a* as unit vectors in the directions of increasing Saturnicentric distance 
r, colatitude 0, and east longitude ty, the magnetic field vector if for the present model 
includes only terms proportional to the five coefficients in Table 1. The results for the field 
components are 


— G,° (c*6e) - G, 

s - V- 3” *■ [ 3 ( +33 [i 

+ - G-, 


( 11 ) 


( 12 ) 


and • (13) 

Note that the external terms (proportional to Gf and G/ ) correspond to a field of magnitude 
[(G i ) 2 +<G i 1 ) 2 ] 1/2 = 1 .005x1 CT 4 G, independent of distance r. In regions in which the magnetic 
lines of force do not connect to this external field, energetic charged particles may be 
trapped, in which case their motion can be approximately separated into three components: 

(1 ) circulation about the field line, 

(2) bounce along the field line, and 

(3) drift around the planet. 


These motions are listed in order of increasing path length and period (decreasing 
frequency), and are in addition to those of the field with Saturn’s rotation and orbit. For an 
individual particle of rest mass m 0 the constants associated with the circulation are the 
kinetic energy E and the magnetic moment M (sample units J/T or MeV/G), as related by 






( 14 ) 
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Here the momentum p and mirror point magnetic field strength B m are given by 

•f>V-E(e>2-n,c*0 (15) 

and 3^ = B (an**') ^ • (16) 


The constants associated with the bounce motion are B m and the adiabatic invariant / given 

by 

r*y»Z '/i 

. (17) 

This integral over path length proceeds along the field line between the mirror points, at each 
of which B=B m , a=90 deg, and the momentum f is perpendicular to i? m . Elsewhere the 
pitch angle a is in the range (0, 90 deg), having a minimum value a e given by 

(W 4 e ) * (18) 


at the magnetic equator (minimum field strength B e between s mi and s m2 along the field line). 
The constant associated with the drift motion is the magnetic shell parameter L , related to / 
and B e as discussed in Section 1 (para. 2). 

For collisionless steady state distributions, the phase space density g (sample SI units 
(J s) -3 ) for a given species of identical particles is constant along the trajectory of every 
individual particle. For present purposes we approximate this relationship using instead 
equal differential intensity (cf. eq. 2) 


i=P 2 g (19) 

at corresponding pitch angles (cf. eq. 16 and 18) everywhere along those field lines 
traversed by the guiding center of the particle circulation. This approximation allows us to 
describe the distribution everywhere in the stable trapping region by specifying only the 
distribution in energy E and pitch angle a e at the magnetic equator and at each value of L , 
as in paragraph 1 of Section 3. 
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APPENDIX B. SPECIES IDENTIFICATION AND OTHER DETECTOR ISSUES 

For some detectors the responses and profiles observed in Saturn’s radiation belts imply 
that the particle species being recorded differs from that anticipated based on experience in 
the Jupiter flybys. The literature includes several discussions of species identification issues 
(see following paragraphs). Many profiles which might otherwise have been used in the 
models have been ignored because these discussions bring their species identification into 
question, or suggest that it’s variable. For the present models, all profiles attributed to ions 
in the original literature are used as if due to protons alone; note that MD80 (p. 5827) and 
SI80 (p. 5738 and p. 5744, Fig. 16) claim that protons are the dominant ion species in 
diverse regions of Saturn’s radiation belts (for a contrasting view, see KR83, pp. 8884 to 
8885). 

Appendix B (pp. 5692 to 5693) of VA80b attributes the GTT countrates A and B 
exclusively to electrons, but C, D, and G to protons (unlike the situation at Jupiter), as 
reflected in Tables 3 and 4. Exceptions occur for troughs in the C and D countrates, and for 
the peaks in the G countrate inside L =4, so these latter features are not used in the model. 
Confirmation for the conversions of the proton rates corresponding to thresholds £>80 MeV 
(first 11 rows of Table 4) is provided by the fact that near L= 2.7 an equatorial flux 
7=2.6xio 4 cm -2 s _1 (VA84, p. 286) yields the C and D peaks well, whereas J=2.0xi0 4 cm _2 s _1 
reproduces the TRD detector M peaks (with the conversions in Table 4). The TRD detector 
M profile ratios confirm that their rates respond to the same species as for GTT C and D 
(protons in the peaks, not in the troughs), as discussed in FM80 (pp. 5804 to 5805). One 
profile from MD80 (p. 5825, Fig. 21, curve C) resembles those just discussed but was not 
used in the model because no absolute conversion is cited and the response is to a mixture 
of proton energies. 

Similarly P11 CPT profiles (for electrons, 7 to 17 MeV) from SI80a (pp. 413 and 414, Fig. 
3 and 5B) and from SI80b (p. 5742, Fig. 13) have not been included in the electron model, 
considering challenges in VA80c (p. 5716) and C089 (p. 5-C-9, para. 3). This last reference 
contains numerous other discussions of particle and energy identifications which have mostly 
been heeded in the present models. An exception is made in that the model does use P1 1 
CPT prbton data (0.5 to 1.8 MeV, several rows near the bottom of Table 4) which are also 
challenged by C089 (p. 5-C-10, bottom), with the result that these data are nor well matched 
(see residuals in Table 8). The model also uses V02 LECP electron data from VA84 (p. 307, 
Fig. 20) in spite of a challenge to its validity (H089, p. 5-D-12, last para.) inside 4 R s ] these 
data consistently exceed the model predictions in this region. 
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Three published spectra were used to help confirm the profile conversions {R/J in Tables 
3 and 4). For electrons, a spectrum at L= 9.56 (VA84, p. 308, Fig. 21) was found to be 
consistent with the conversions used for countrates read from four V02 LECP profiles listed 
in Table 3 (from VA84, p. 307, Fig. 20). For protons, a spectrum at 7.81 R s (MD80, p. 5816, 
Fig. 4) was found to be consistent with the conversion used for differential intensities from 
three P1 1 CRT profiles listed in Table 4 (from MD80, p. 5815, Fig. 3, and from TR80, p. 422, 
Fig. 1). However, also for protons, a spectrum at L=7.15 (VA84, p. 297, Fig. 11) did not 
seem to be consistent with conversions intended for countrates from two V01 LECP profiles 
(KR83, p. 8877, Fig. 6, curve 2 from top, and HA83, p. 8908, Fig. 3a, top curve), so these 
two profiles were not used in the model. 

Adjustments to the scale of the published profiles were made in two cases, both for 
electrons, one of them described in note [j] to Table 3. The other concerns the P11 TRD 
CDC data (Table 3, row 4), for which the original current scale (FM80, p. 5804, Fig. 1) 
appears to increase downwards; the intensity profile (FI80, p. 427, Fig. 2A) and the 
geometric factor 7.3xlCT 14 A cm 2 s sr (FM80, p. 5805, Table 1) suggest that the 10~ 7 A tic 
mark on the current scale is correct (the others should represent 10 -9 and l(T n A). The 
conversion adopted in Table 3 (/?//=6.3xio -15 A cm 2 s) is the ratio 

(4jc srr^lCT 7 - 82 AK1.9X10 5 cm -2 s _1 sr _1 r', where the quantities in brackets are the peak 
inbound current (FM80, p. 5804, Fig. 1) and the peak inbound integral intensity (FI80, p. 427, 
Fig. 2A). 
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APPENDIX C. PROCESSES OF MODEL DEVELOPMENT 

The following paragraphs detail various aspects of the model development, particularly 
including the selection, measurement, and filtration of the data profiles (cf. Sec. 2 and Tables 
3 and 4) and the iteration of the model parameter values (cf. Sec. 4, para. 3 and 4). 

To choose from among the many data profiles considered for possible inclusion in the 
model, a selection process was used based roughly on the following considerations: 

(1) unpublished data were not sought, on the assumption that the published data may be 
the most secure, in terms of careful reduction and analysis; 

(2) figures which display the spacecraft data as spectra or as phase space densities were 
not used, on the assumption that most values therein are already represented in the 
profiles of countrate, flux, or intensity (but spectra were used to confirm conversions in 
some cases); 

(3) data were used only when an unambiguous, constant conversion between the profile 
quantity and flux (in absolute units, here cm _2 s _1 ) could be determined from 
accompanying, published information (see eq. 7); 

(4) consistent with item (3), very noisy data, and data requiring variable countrate 
corrections (e.g., for dead-time), were mostly avoided; and 

(5) the same, identifiable species (e.g., electrons or protons), at constant threshold energy 
(again see eq. 7), is responsible for that portion of the profile used in the model. 

With respect to consideration (5), a brief discussion of species identification within the data 
set appears in Appendix B. As a result of these ideas, a total of 239 profiles identified in the 
literature for potential use in the models was narrowed down to the 31 cited for electrons in 
Table 3 and the 32 cited for protons in Table 4. 

In addition to the use of countrate, flux, or intensity (detector current in one case) as the 
profile quantity, the several sources display data variously as functions of Earth-received 
time, spacecraft time, distance from Saturn, or magnetic shell parameter L. In most cases, 
profiles along the inbound and outbound trajectory legs are distinctly displayed. For some 
detectors or channels, more than one published profile applies (e.g., rows 3 through 5 in 
Table 4), in which case the model uses profile points only from that source having the 
greatest resolution in time (or distance or l) available at each trajectory point (see note [a ] 
to Tables 3 and 4). Unit weight was assigned to each detector or channel, with the very few 
exceptions cited h the notes to Tables 3 and 4. 
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Positions ( X and Y in cm, using as origin the lower left corner of the figure frame) were 
measured for enough points along each profile to capture all the maxima, minima, gaps, and 
other important features (plateaus, shoulders, etc.). These positions were then converted 
linearly to the quantities with which the figure scales were labelled (commonly log R for the 
profile quantity), and the resulting digitized data were stored in computer files. The data from 
these files were replotted (by machine) on the same scales as the originals to enable visual 
verification that the original profiles were satisfactorily represented by the digitized data files. 
These files were then compared with others describing the three spacecraft trajectories (as 
generated from the elements in Table 2) to extract those values of L at which the profiles 
have maxima, minima, and end-points. Coincidence in L of the maxima and minima from 
the diverse profiles was generally satisfactory, enabling a selection of suitable values of L for 
use in the models (column 1 of Tables 5 and 6). Then, at each such value of L, using the 
times and distances in the spacecraft trajectory files, values were interpolated from the 
profile data files yielding a data point log R D (e.g., countrate, flux or intensity) for each of 
several detectors (identified by local magnetic field B , energies E a and E b , and conversion 
RIJ). 

The iteration for model parameter values was completed independently at each value of 
L , using the considerations in Section 4 (para. 3 and 4). Additional features incorporated in 
the interactive computer program used for the iteration included 

(1 ) control over which parameters to iterate in a given attempt, and within what range of 
values, as well as allowing any b n value to be slaved to a n when desirable (see eq. 3); 

(2) capability to start with a pure power law (by taking E 2 and £ 3 respectively much smaller 
and much larger than the range of energies E a and E b in the data set) and to retain the 
values so obtained for a 2 and b 2 (see Sec. 3, para. 2, item 1); 

(3) occasional, temporary reweighting of the data to ensure that small residuals were 
maintained for some key countrates during parameter variation (however, the level 
weighting cited in the text was always reintroduced at the end); and 

(4) progressive reduction of A, from some large value (in which case eq. 9 and 1 0 yield Y 
as a standard RMS residual) to 0.3 (see end of para, following eq. 10), such that the 
"hopeless” data points are successively deemphasized in the iteration (allowing 
concentration on those having |A| near A 0 , as being most amenable to further residual 
reduction). 

With regard to feature (4) above (cf. eq. 8 through 10), note that A=» leads only to Z=2A} or 
to Z 1/2 =0.424 when A o =0.3; cf. the value T=0.35 for protons at L =2.52 (Fig. 3, last column), 
which represents by far the least satisfactory result in the model iteration process. 
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Table 1 . Parameter values for Saturn used in the model calculations. 

QUANTITY 

VALUE 

SOURCE 

Radius 

Gravitation 

Rotation 

Magnetic field coefficients 

i?s=6.033xl0 7 m 
GM S =3 .794x10 16 m 3 /s 2 
co=810.76 deg/day 
g° =0.21 535 G ^ 

g2 =0.01 642 G 
#3=0.02743 G 
Gf =-0.0001 G 

Gi =-0.00001 G J 

C084, p. 354 
CE89 
C084, p. 360 

| C084.P. 363. Table II 
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Table 2. Circumstances of the flybys for three spacecraft at Saturn. The Pioneer 1 1 data 
are from various sources in the Science issue which includes VA80a (and note 6 
on p. 421 includes the one-way-light-time 86.34 min, Satum-to-Earth, used for 
profile data interpretation). 


Spacecraft 


Pioneer 1 1 

Voyager 1 

Voyager 2 

Perisaturn 

date 

1979 Sept. 1 

1980 Nov. 12 

1981 Aug. 26 


day of year 

244 

318 

238 


spacecraft time (UT) 

16:31:00 

23:45:43 

03:24:05 


distance, r l {R s ) 

1.3497 

3.0522 

2.67075 


argument, oox (deg) 

-88.7 

-75.7 

153.2 


west longitude, j2, (deg) 

162.1 

179.9 

322.9 

Ascending node 

west longitude, ^ (deg) 

139 

265 

150 

Eccentricity 

e 

1.1148 

2.1078 

1.48587 

Inclination 

i (deg) 

6.55 

41.17 

29.51 

Source 


See caption 

CE80, CE89 

W081 , CE89 
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Table 3. Identification of data profiles used for the Saturn electron model. Detailed notes 
appear on the following page. 


SPACECRAFT, 

SUBSYSTEM 

[a] 

REFERENCE, 

PAGE 

FIGURE, 

LINE 

[ft] 

CHANNEL 

ENERGIES 

Ea>Eb 

(MeV) 

[c] 

CONVERSION 

R/J 

(cm 2 ) 

[d] 

BACK- 

GROUND 

\ogR B 

[e] 

NOTES 

V02 LECP 

KR82, 573 

2 C, 2 

EAB 10 

10 , 20 

1.03(-3) 

• 2.0 

If] 

P 11 CPT 

SI80b, 5749 

21 c, 1 

ECD 

3.4 

1.00 

4.2 

[£] 


SI80b, 5742 

13.4 

ECD 

3.4 

1.00 

4.2 

[£] 

P 1 1 TRD 

FM80, 5804 

lab, 4 

CDC 

2.0 

6.3(-15) 

0.0 

[ft] 

P 1 1 CPT 

SI80b, 5740 

11,2 

MT 

2.0, 7.0 

3. 42 (-2) 

- 1.6 

[i] 

V02 LECP 

KR82, 573 

2C, 1 

ESA 0 

1.5 

9.55(-4) 

0.0 

U] 

P 1 1 CRT 

MD80, 5824 

18,4 


1 . 1 , 2.0 

0.088 

1.0 

[ft] 


TR80, 422 

2,3 


1 . 1 , 2.0 

0.088 

1.0 

[ft] 

V02 LECP 

VA84, 307 

20,5 

E 7 O 8 -E 7 O 9 

0.853, 1.2 

1.35(-5) 

-1.5 

[/] 

P11 CRT 

MD80, 5824 

18,3 


0 . 8 , 1.1 

0.265 

2.0 

[ft] 


TR80, 422 

2,2 


0 . 8 , 1.1 

0.265 

2.0 

[ft] 

P 11 GTT 

VA80a, 419 

5,2 

B 

0.56, 21 

1.27(-3) 

-0.4 

[m] 


VA80c, 5710 

1,2 

B 

0.56, 21 

1.27(-3) 

-0.4 

[m] 


VA84, 300 

14,2 

B 

0.56, 21 

0.08 

2.0 

In] 

V02 LECP 

VA84, 307 

20,4 

E'j07-E'y08 

0.48, 0.853 

2.79(-4) 

1.0 

[l] 

P 11 TRD 

FI80, 429 

5,1 

El 

0.45 

1.64(-3) 

0.2 

[0] 


FM80, 5804 

lab, 3 

E3 

0.46 

4.54(-4) 

0.0 

[ P ] 

P 1 1 CRT 

MD80, 5824 

18,2 


0.43, 0.8 

0.215 

2.0 

[ft] 

V02 LECP 

VA84, 307 

20,3 

E 7 O 6 -E 7 O 7 

0.252, 0.48 

6.45(-4) 

1.7 

[/] 


AR83, 8899 

12,1 

E 7 O 6 -E 7 O 7 

0.252, 0.48 

6.4S(-4) 

0.5 

[f] 

V 01 LECP 

KR81 , 226 

1 A, 2 

E(}05 

0.2, 0.5 

1.59(-4) 

- 1.0 

If] 

P 1 1 CRT 

MD80, 5824 

18. 1 


0.16, 0.43 

0.295 

3.0 

[ft] 


TR80, 422 

2 , 1 


0.1,0.43 

0.241 

3.0 

[ft] 

P 11 TRD 

FI80, 427 

2C, 2 

E1-1.25E2 

0.16, 0.25 

0.856 

4.0 

[?] 

V 01 LECP 

CH85, 508 

2, 3&4 

EP04 

0.14, 0.2 

4.77(-4) 

-0.5 

If] 


KR83, 8877 

6,5 

Ep04 

0.13, 0.2 

4.77(-4) 

- 1 .0 

[f] 

V02 LECP 

VA84, 307 

20,2 

EP04 

0.112, 0.183 

3.1 (-4) 

1.0 

[l] 


AR83, 8899 

11 * 2 

E(304 

0.112, 0.183 

3.1 (-4) 

1.0 

[f] 

P 11 GTT 

VA80a, 419 

5, 1 

A 

0.04, 21 

1.37{-3) 

-0.3 

[m] 


VA80c, 5710 

1 , 1 

A 

0.04, 21 

1.37(-3) 

-0.3 

[m] 


VA84, 300 

13, 1&2 

A 

0.04, 21 

0.08 

2.7 

[n] 
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NOTES TO TABLE 3: 

[a] The profiles appear in order of decreasing threshold energy E a (see column 5, note [c]); a profile lacking 
an entry in column 1 is used in the model only for trajectory points at which no value for countrate R (note 
[, d ]) is available from the profile immediately above. 

[b ] Second entry in this column is the curve placement, counting down from the top of the figure cited. 

[c] Where two entries appear, the response (next column) is assumed constant between them, and zero 
elsewhere; a single entry represents the smaller (threshold E a ) of these two energies, and the other (E b ) is 
assumed infinite. 

[, d ] Units cm -2 assumes R is a countrate (units s" 1 ); exceptions are cited in the notes, and the entry in this 
column is always the ratio of the profile quantity R (countrate, intensity, flux, etc.) to the local omni flux J 
(see text following eq. 7). Parentheses enclose powers of ten. 

[e ] This entry is estimated from the noise or limit of the profile in the figure, in the same units as R in the 
preceding column (see note [d]). 

[f ] Energies from profile label, conversion is eG in cm 2 sr (KR83, p. 8874, Table 1 ) divided by 4 jc sr. 

[g ] Energy from figure label, profile is flux J , no conversion needed. 

[h ] Energy from source (FM80, p. 5805, Table 1), profile is current, conversion is as discussed in Appendix B 
(last para.). 

[i ] Energies from profile label, conversion is 0.43 cm 2 sr (SI80b, p. 5757, Table A1) divided by 47c sr. 

[y ] Energy from profile label, conversion is eG-0.12 cm 2 sr (KR83, p. 8874, Table 1) divided by 4rc sr and by an 
additional factor 10 to secure agreement with a corresponding profile (KA82, p. 1144, Fig. 2, curve 3 from 
top) and with Pioneer 1 1 detectors at nearby energies. 

[k ] Energies from figure caption (or label), profile is differential intensity i , conversion is (4rt sr) -1 (Ej, — E a ) . 

[/ ] Same as note [f ], but conversion also confirmed by comparison with spectrum at L -9.56, as discussed in 
Appendix B (para. 4). 

[m ] Energies from source (VA80a, p. 416, Table 1), conversion is Q from same table. 

[n ] Energies as in note [m ], profile is integral intensity /, conversion is (4 tc sr) -1 . 

[o ] Energy from profile label, comversion is (/? U )j(R i IR 3), where subscript 1 is for channel El , subscript 3 is 
for channel E3, R \ and R 3 are the peak countrates at L =2.7, and conversion (R //) 3 is derived in note [p ]. 

\p ] Energy from source (FM80, p. 5805, Table 1), conversion is 0.0057 cm 2 sr (same table) divided by 4n sr. 

[q ] Energies from source (FI80, p. 426, Table 1), profile is differential intensity i , conversion is 
(47C sr)-\E b -E a y\ 
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Table 4. Identification of data profiles used for the Saturn proton model. Detailed notes 
appear on the following page. 


SPACECRAFT, REFERENCE, 

SUBSYSTEM PAGE 

[a] 

FIGURE, 

LINE 

[b] 

CHANNEL 

ENERGIES 

EaEb 

(MeV) 

[c] 

CONVERSION 

R/J 

(cm 2 ) 

[d] 

BACK- 

GROUND 

log/? s 

[e] 

NOTES 

P1 1 TRD 

FI80, 429 

5,3 

Ml 

80 

0.038 

0.0 

if) 


FM80, 5804 

lab, 5 

M3 

80 

0.012 

-1.0 

[g] 

P11 GTT 

VA80a, 418 

4, 2&4 

D 

80 

0.043 

-0.5 

[h] 


VA80c, 5710 

1.4 

D 

80 

0.043 

-0.5 

[i] 


VA80a, 417 

3, 2&4 

D 

80 

0.043 

-0.5 

[h] 

P11 GTT 

VA80a, 419 

5,3 

C 

80 

0.122 

-0.3 

[A] 


VA80a, 418 

4, 1&3 

C 

80 

0.122 

0.0 

[h] 


VA83, 6913 

4,1 

C 

80 

0.122 

0.0 

111 


VA83, 6913 

5,1 

C 

80 

0.122 

0.0 

V) 


VA80c, 5710 

1.3 

C 

80 

0.122 

-0.3 

[i] 


VA80a, 417 

3, 1&3 

C 

80 

0.122 

0.0 

[h] 

V02 CRS 

SM83, 8925 

3.2 


63, 160 

8.2(-4) 

-3.0 

[*] 

V02 LECP 

KR82, 573 

2C, 3 

PAB12 

54, 87 

0.004 

-2.0 

[/] 

V02 CRS 

SM83, 8925 

3, 1 


48, 63 

5.31 (-3) 

-2.5 

[k] 


SM83, 8924 

2,1 


48, 63 

1.85 

0.0 

[m] 

P1 1 CPT 

SI80b, 5742 

14,2 

FCD 

35 

2.5(-5) 

-1.5 

[«] 

V02 CRS 

V082, 579 

3, 1 

A 

27, 160 

0.08 

-0.8 

[o] 

P1 1 CRT 

MD80, 5815 

3,4 


2.1, 3.1 

0.08 

-2.0 

Ik) 

P11 CRT 

MD80, 5824 

19, 1 


1.56, 5.1 

0.08 

-2.0 

[0} 

P1 1 CRT 

TR80, 422 

1,3 


1.1,2.15 

0.076 

0.0 


P11 CRT 

TR80, 422 

1.2 


0.75, 1.2 

0.177 

1.0 

IP] 

P11 GTT 

VA80a, 417 

2,1 

G 

0.61,3.41 

3.5(-3) 

-0.8 

[h] 


VA84, 290 

5, 1&2 

G 

0.61 , 3.41 

0.08 

0.2 

[o] 

V02 LECP 

VA84, 296 

10,2 

PL06 

0.54, 0.99 

0.009 

-3.0 

[D 

P1 1 CRT 

MO80, 5815 

3,2 


0.53, 0.74 

0.379 

1.0 

[Pi 

P11 CPT 

SI80b, 5748 

20A2, 1 

LI L2 

0.5, 1.8 

0.039 

0.3 

lq] 


SI80b, 5748 

20C2, 1 

LI [2 

0.5, 1 .8 

0.039 

0.3 

lq] 


SI80b, 5742 

14,1 

L1L2 

0.5, 1.8 

0.039 

0.3 

lq] 


SI80b, 5743 

15, 1&2 

L1L2 

0.5, 1.8 

0.039 

0.3 

lq 1 

P11 CRT* 

TR80, 422 

1.1 


0.2, 0.75 

0.145 

2.0 

[k] 

V01 LECP 

CH85, 508 

1.3&4 

PL04 

0.139, 0.22 

3.2(-3) 

-0.5 

[D 

V02 LECP 

AR83, 8895 

3,2 

PL04 

0.137,0.215 

0.009 

2.0 

ID 
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NOTES TO TABLE 4: 

[a ] The profiles appear in order of decreasing threshold energy E a (see column 5, note [c ]); a profile lacking 
an entry in column 1 is used in the model only for trajectory points at which no value for countrate R (note 
[rf]) is available from the profile immediately above. 

[b ] Second entry in this column is the curve placement, counting down from the top of the figure cited. 

[c] Where two entries appear, the response (next column) is assumed constant between them, and zero 
elsewhere; a single entry represents the smaller (threshold E a ) of these two energies, and the other (£j, ) is 
assumed infinite. 

[4] Units cm -2 assumes R is a countrate (units s _1 ); exceptions are cited in the notes, and the entry in this 
column is always the ratio of the profile quantity R (countrate, intensity, flux, etc.) to the local omni flux J 
(see text following eq. 7). Parentheses enclose powers of ten. 

[e ] This entry is estimated from the noise or limit of the profile in the figure, in the same units as R in the 
preceding column (see note [4]). 

| [f ] Energy from profile label, conversion is geometric factor in cm 2 (FI80, p. 426, Table 1 ). 

[g ] Energy from source (FM80, p. 5805, Table 1), conversion is geometric factor in cm 2 (same table; note 
disagreement with FI80, p. 426, Table 1). 

[h ] Energy from source (VA80a, p. 416, Table 1), conversion is Q from same table. 

[i ] Energy from source (VA80b, p. 5682, Table 2), conversion is Q from same table. 

[7 ] Energy as in note [h ], conversion from source (VA83, p. 6912, eq. 1). 

[k ] Energies from figure caption (or label), profile is differential intensity i , conversion is (451 sr) -1 (£ fe —E a ) . 

[/ ] Energies from profile label, conversion is zG in cm 2 sr (KR83, p. 8874, Table 1 ) divided by 4tt sr. 

[m] Energies as for preceding profile, conversion is (4ic sr) -1 (Ej, —E a ) _1 (/? 3/1 3) where log/? 3-3.75 and 
logi 3-1.21 are peak rates from this and the preceding profile. 

[n ] Energy from profile label, conversion from source (SI80b, p. 5757, Table A1), weight is 2.0 to compensate 
for indistinguishability of the inbound and outbound profiles. 

[o ] Energies from figure caption, profile is integral intensity / , conversion is (47C sr) -1 . 

[ p ] Same as note [k ], but conversion also confirmed by comparison with spectrum at 7.81 R s , as discussed in 
Appendix B (para. 4). 

[q ] Energies from profile label, conversion is geometric factor in cmV (SI80b, p. 5757, Table A1 ) divided by 4ic 
sr. 
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Table 5. Parameter values for the Saturn electron model . See model description 
accompanying equations (1) through (5). 


L 

a 0 


a 2 

a 3 

K 

b i 

b 2 

t>3 

e 2 

E 3 

2.30 

4.74 

0.00 

0.00 

6.47 

4.74 

0.00 

0.00 

6.47 

0.60 

1.50 

2.33 

5.02 

0.00 

0.00 

4.01 

5.02 

0.00 

0.00 

4.01 

0.60 

1.50 

2.34 

4.91 

0.00 

0.00 

4.86 

4.91 

0.00 

0.00 

4.86 

0.60 

1.50 

2.36 

4.98 

0.00 

0.00 

4.41 

4.98 

0.00 

0.00 

4.41 

0.60 

1.50 

2.38 

5.00 

0.00 

0.00 

5.37 

5.00 

0.00 

0.00 

5.37 

0.60 

1.50 

2.46 

5.49 

0.00 

0.23 

3.00 

5.49 

0.00 

0.23 

3.00 

0.60 

2.00 

2.56 

5.56 

0.00 

0.34 

3.00 

5.56 

0.00 

0.34 

3.00 

0.60 

2.00 

2.70 

5.51 

0.00 

0.27 

4.48 

5.51 

0.00 

0.27 

4.48 

0.60 

2.00 

2.83 

5.52 

0.00 

0.25 

3.11 

5.52 

0.00 

0.25 

3.11 

0.20 

2.00 

2.93 

5.46 

0.00 

0.28 

3.42 

5.46 

0.00 

0.28 

3.42 

0.20 

2.00 

3.04 

5.38 

0.00 

0.27 

4.96 

5.38 

0.00 

0.27 

4.96 

0.20 

2.00 

3.08 

5.31 

0.00 

0.29 

5.23 

5.48 

0.00 

0.29 

5.23 

0.20 

2.00 

3.20 

5.03 

0.00 

0.20 

5.48 

5.74 

0.00 

0.20 

5.48 

0.20 

2.00 

3.40 

5.50 

0.00 

0.00 

4.27 

4.90 

0.00 

0.61 

4.27 

0.20 

2.00 

3.90 

4.88 

0.00 

0.20 

5.71 

4.72 

0.21 

0.78 

5.71 

0.20 

2.00 

4.30 

4.67 

0.00 

0.00 

4.29 

4.50 

0.00 

0.56 

4.29 

0.10 

1.00 

4.44 

4.74 

0.00 

0.00 

4.23 

4.12 

0.00 

0.68 

4.23 

0.10 

1.00 

4.96 

3.58 

0.00 

0.34 

5.74 

5.30 

0.00 

0.11 

5.74 

0.10 

1.00 

5.41 

3.50 

0.00 

0.27 

5.26 

5.38 

0.00 

0.11 

5.26 

0.10 

1.00 

5.87 

4.57 

0.00 

0.05 

6.67 

3.73 

1.01 

0.68 

6.67 

0.10 

1.00 

6.20 

3.26 

0.00 

0.00 

4.22 

4.60 

0.26 

0.00 

4.22 

0.08 

0.50 

6.50 

3.95 

0.00 

0.06 

3.00 

4.83 

0.00 

0.19 

3.00 

0.08 

0.50 

7.20 

3.33 

0.00 

0.00 

6.10 

3.20 

0.21 

0.66 

6.10 

0.08 

0.50 

8.00 

3.42 

0.00 

0.00 

5.50 

2.56 

0.16 

2.45 

5.50 

0.08 

0.50 

9.50 

3.02 

0.00 

0.83 

6.00 

2.11 

0.00 

2.22 

6.00 

0.08 

0.50 

11.00 

3.64 

0.00 

1.17 

3.01 

1.98 

0.00 

3.78 

3.01 

0.08 

0.50 

13^.00 

2.09 

2.23 

0.00 

7.59 

2.09 

2.23 

0.00 

7.59 

0.08 

0.50 
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Table 6. Parameter values for the Saturn proton model. See model description 
accompanying equations (1) through (5). 


L 

a» 

a i 

a 2 

a 3 

b 0 

b i 

b 2 

b 3 

E 2 

E 3 

2.30 

1.86 

0.00 

0.00 

4.57 

1.86 

0.00 

0.00 

4.57 

1.00 

60.0 

2.32 

2.07 

0.00 

0.00 

4.53 

2.07 

0.00 

0.00 

4.53 

1.00 

60.0 

2.34 

2.17 

0.00 

0.00 

6.09 

2.17 

0.00 

0.00 

6.09 

1.00 

60.0 

2.43 

3.18 

0.05 

0.00 

3.05 

3.18 

0.05 

0.00 

3.05 

1.00 

60.0 

2.52 

2.14 

1.85 

0.00 

3.00 

2.14 

1.85 

0.00 

3.00 

1.00 

60.0 

2.67 

3.80 

0.75 

0.00 

3.00 

3.80 

0.75 

0.00 

3.00 

1.00 

60.0 

2.77 

3.37 

1.28 

0.00 

3.00 

3.37 

1.28 

0.00 

3.00 

0.50 

60.0 

2.82 

3.26 

2.63 

0.00 

3.00 

3.26 

2.63 

0.00 

3.00 

0.30 

60.0 

2.86 

3.14 

2.93 

0.00 

3.00 

3.14 

2.93 

0.00 

3.00 

0.30 

60.0 

3.10 

0.37 

0.00 

0.00 

3.00 

0.67 

4.87 

6.27 

3.00 

1.00 

60.0 

3.20 

2.91 

0.29 

0.00 

3.00 

-0.25 

6.63 

0.00 

3.00 

1.00 

60.0 

3.31 

2.86 

0.02 

0.02 

3.00 

0.29 

8.30 

0.00 

3.00 

0.50 

60.0 

3.44 

2.80 

0.98 

0.00 

3.00 

0.46 

7.41 

0.00 

3.00 

0.50 

60.0 

3.71 

2.16 

1.77 

0.00 

3.00 

0.40 

7.10 

0.00 

3.00 

0.50 

60.0 

3.99 

1 

O 

—A. 

CO 

0.00 

0.26 

3.00 

0.58 

7.70 

0.52 

3.00 

0.50 

60.0 

4.13 

1.29 

0.00 

0.00 

3.00 

-0.10 

9.00 

0.15 

3.00 

0.50 

60.0 

4.27 

0.63 

0.46 

0.29 

3.22 

0.16 

8.55 

0.00 

3.22 

0.50 

60.0 

4.38 

0.61 

5.05 

0.20 

3.12 

-0.03 

6.93 

0.00 

3.12 

0.50 

60.0 

4.49 

0.47 

5.25 

0.21 

3.55 

-0.06 

7.15 

0.00 

3.55 

0.50 

60.0 

5.01 

0.68 

1.45 

0.00 

4.41 

-1.13 

6.73 

6.85 

4.41 

0.50 

0.80 

5.36 

1.77 

0.00 

0.00 

3.00 

-0.98 

6.78 

8.25 

3.00 

0.50 

0.80 

6.20 

2.32 

3.32 

0.00 

3.08 

0.10 

5.51 

2.99 

3.08 

0.50 

0.80 

6.45 

2.40 

3.34 

0.00 

3.00 

0.06 

5.81 

2.66 

3.00 

0.50 

0.80 

6.70 

2.46 

2.59 

0.52 

4.72 

0.30 

5.83 

0.12 

4.72 

0.50 

0.80 

7.63 

2.46 

3.55 

0.25 

3.85 

0.81 

4.90 

0.61 

3.85 

0.50 

0.80 

8.40 

2.31 

6.28 

1.73 

3.35 

1.07 

0.48 

1.73 

3.35 

0.30 

0.80 

10.30 

1.90 

6.52 

1.81 

4.11 

0.74 

0.00 

0.85 

4.11 

0.30 

0.80 

12.00 

1.24 

0.00 

6.35 

3.94 

1.24 

0.00 

6.35 

3.94 

0.30 

0.80 

13.00 

1.22 

0.00 

6.72 

3.31 

1.22 

0.00 

6.72 

3.31 

0.30 

0.80 
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Table 7. Summary of range in L and residuals (see eq. 8) for the Saturn electron model. 

See Table 3 (column 1 and note [a]) for association with the data profiles by 
detector. 


SPACECRAFT, 

MINIMUM 

MAXIMUM 

MINIMUM 

MEAN 

MAXIMUM 

RMS 

SUBSYSTEM 

L 

L 

A 

A 

A 

A 

V02 LECP 

2.83 

5.41 

-0.16 

0.01 

0.14 

0.09 

P11 CPT 

2.34 

2.38 

-0.21 

-0.03 

0.19 

0.14 


2.33 

5.41 

-0.07 

0.72 

1.23 

0.81 

P11 TRD 

2.36 

6.20 

-0.66 

-0.09 

0.17 

0.20 

P11 CPT 

9.5 

13.0 

-3.23 

-1.76 

-0.57 

2.04 

V02 LECP 

2.83 

5.41 

-1.71 

-0.13 

0.18 

0.46 

P11 CRT 

2.30 

11.0 

-0.81 

0.18 

1.05 

0.41 


13.0 

13.0 

-0.02 

0.35 

0.71 

0.50 

V02 LECP 

2.83 

9.50 

-1.35 

0.50 

1.42 

0.93 

P11 CRT 

2.30 

11.0 

-0.48 

0.09 

0.90 

0.26 


13.0 

13.0 

-0.03 

0.39 

0.80 

0.57 

P11 GTT 

2.30 

2.46 

-0.26 

0.10 

0.52 

0.28 


2.70 

3.40 

-0.04 

0.03 

0.13 

0.06 


2.56 

13.0 

-0.27 

0.02 

0.21 

0.13 

V02 LECP 

2.83 

9.50 

0.02 

0.41 

0.84 

0.45 

P1 1 TRD 

2.33 

3.20 

-1.35 

-0.22 

0.09 

0.46 


2.70 

6.20 

-0.18 

-0.05 

0.16 

0.11 

P1 1 CRT 

2.30 

11.0 

-0.56 

-0.12 

0.28 

0.21 

V02 LECP 

2.93 

9.50 

-0.34 

0.36 

0.87 

0.45 


2.83 

11.0 

-0.20 

0.07 

0.62 

0.19 

V01 LECP 

4.44 

13.0 

-1.63 

-0.05 

0.95 

0.60 

P11 CRT 

4.96 

11.0 

-0.21 

0.07 

0.35 

0.18 


13.0 

13.0 

-0.15 

-0.04 

0.08 

0.12 

P1 1 TRD 

5.41 

13.0 

-0.69 

-0.17 

0.18 

0.28 

V01 LECP 

11.0 

13.0 

0.08 

0.46 

0.92 

0.55 


4.44 

9.50 

-0.11 

0.36 

1.24 

0.58 

V02 LECP 

2.83 

9.5 

-0.07 

0.27 

0.55 

0.34 


2.83 

7.20 

-0.38 

-0.09 

0.05 

0.14 

P11 GTT 

2.33 

2.46 

-0.07 

0.14 

0.43 

0.23 


2.70 

3.40 

-0.20 

-0.13 

-0.05 

0.14 


2.36 

13.0 

-0.46 

-0.11 

0.05 

0.17 
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Table 8. Summary of range in L and residuals (see eq. 8) for the Saturn proton model. 

See Table 4 (column 1 and note [a]) for association with the data profiles by 
detector. 


SPACECRAFT, 

MINIMUM 

MAXIMUM 

MINIMUM 

MEAN 

MAXIMUM 

RMS 

SUBSYSTEM 

L 

L 

A 

A 

A 

A 

P1 1 TRD 

2.43 

3.31 

-0.93 

-0.20 

0.15 

0.41 


2.77 

4.13 

-1.71 

-0.07 

0.76 

0.59 

P11 GTT 

2.43 

2.82 

-0.04 

0.15 

0.73 

0.27 


2.86 

3.44 

-0.12 

0.57 

2.46 

1.11 


2.30 

4.27 

-1.20 

0.17 

1.75 

0.76 

P11 GTT 

2.30 

2.43 

-0.61 

0.00 

0.37 

0.37 


2.52 

2.82 

0.11 

0.25 

0.85 

0.35 


2.82 

3.10 

0.20 

0.98 

2.50 

1.46 


2.86 

3.10 

0.21 

1.40 

2.58 

1.83 


3.20 

3.44 

-0.04 

0.13 

0.23 

0.16 


2.34 

4.27 

-1.06 

0.51 

1.98 

1.01 

V02 CRS 

2.77 

5.36 

-0.82 

-0.15 

0.20 

0.31 

V02 LECP 

2.77 

5.01 

-2.29 

-0.63 

0.06 

0.86 

V02 CRS 

2.77 

5.36 

-1.22 

-0.08 

0.21 

0.32 


2.82 

3.10 

-0.80 

-0.53 

-0.28 

0.57 

P1 1 CPT 

2.67 

3.44 

-0.26 

-0.16 

-0.04 

0.18 

V02 CRS 

4.13 

5.36 

-0.67 

-0.09 

0.08 

0.24 

P1 1 CRT 

3.31 

4.49 

-0.24 

-0.10 

0.03 

0.13 

P1 1 CRT 

2.30 

5.01 

-1.19 

-0.13 

0.13 

0.36 

P1 1 CRT 

2.43 

13.0 

-0.38 

0.27 

1.32 

0.46 

P11 CRT 

5.01 

13.0 

-0.35 

-0.01 

0.38 

0.20 

P11 GTT 

5.36 

7.63 

-0.12 

-0.01 

0.10 

0.08 


3.99 

13.0 

-0.22 

0.10 

0.51 

0.21 

V02 LECP 

3.10 

10.3 

-0.50 

-0.07 

0.19 

0.17 

P11 CRT 

3.31 

5.01 

-0.07 

0.56 

1.17 

0.69 

P11 CPT 

2.43 

2.67 

-1.49 

-0.91 

0.0 

1.12 


2.67 

2.67 

-1.30 

-1.30 

-1.30 

1.30 


3.10 

4.49 

-1.02 

-0.27 

-0.03 

0.44 


4.38 

13.0 

-0.46 

-0.14 

0.24 

0.21 

P1 1 CRT 

2.43 

13.0 

-0.44 

0.12 
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Figure 1. Trajectories of Saturn flyby spacecraft in magnetic coordinates (reproduced from 
Fig. 1c of KR83). The Pioneer 11 trajectory closely follows the curve labelled 
\B eq \ inwards to L=1.43 and 5=0.0801 G, whereas the Voyager 1 and 2 
trajectories follow the dashed and solid curves, respectively. See field model and 
trajectory entries in Tables 1 and 2. 
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Figure 2. Coverage of the data for the Saturn electron model as a function of L. At the left, 
the histogram encloses values of magnetic field ratio B e /B spanned by the data. 
In the center, the histogram encloses values of the threshold energy (E a in Table 
3) spanned by the data (dashed lines are E 2 and £ 3 from Table 5). At the right 
three columns specify the number of data points used, and those having 
residuals limited by the factors 4* 1 and 2 ±1 . The last column displays the modified 
RMS residual at each L (see eq. 10, using A o =0.3). 
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Figure 3. Coverage of the data for the Saturn proton model as a function of L. At the left, 
the histogram encloses values of magnetic field ratio BJB spanned by the data. 
In the center, the histogram encloses values of the threshold energy ( E a in Table 
4) spanned by the data (dashed lines are E 2 and £ 3 from Table 6). At the right 
three columns specify the number of data points used, and those having 
residuals limited by the factors 4 ±1 and 2 ±1 . The last column displays the modified 
RMS residual at each L (see eq. 10, using A„=0.3). 
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Figure 4. Integral intensity for the Saturn electron model at three energy thresholds. As 
evaluated at the magnetic equator, the perpendicular and parallel intensities (l 
and / ft are represented respectively by the solid and dashed lines. Note that the 
curves for £=1 .0 MeV reproduce the values for a Q and b Q in Table 5 (cf. eq. 1 
and 3). The horizontal axis changes scale atL=5. 
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Figure 5. Integral omnidirectional flux for the Saturn electron model at three energy 
thresholds. The fluxes at the magnetic equator (b,/b= 1.0) and at high magnetic 
latitude (B e iB =0) are represented respectively by the solid and dashed lines. The 
horizontal axis changes scale at L= 5. 
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Figure 6. Integral intensity for the Saturn electron model at three sample values of L 
(numbers adjacent to the curves). On the left the dependences on pitch angle a 
are shown for energy thresholds of 0.1, 1.0, and 10 MeV (solid, dashed, and solid 
lines) at each L . On the right spectra are shown for / x (solid lines) and for /„ 
(dashed lines) at each L (/ x and 7 (| are equal for the isotropic distribution at 
L =2.83). 
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Figure 7. Integral omnidirectional flux for the Saturn electron model at three sample values 
of L (numbers adjacent to the curves). On the left the dependences on local 
magnetic field strength B are shown for energy thresholds of 0.1, 1.0, and 10 
MeV (solid, dashed, and solid lines) at each L. On the right spectra are shown 
for the magnetic equator {B e /B= 1.0, solid lines) and for high magnetic latitude 
(B e /B= 0, dashed lines) at each L. 





<*> 

o 


88 


'PROTON IMTE&RAL IMTENSIT^ ^ r 6*rVv 


5217-90-029 


35 


A. R. Hoffman 


Figure 8. Integral intensity for the Saturn proton model at three energy thresholds. As 
evaluated at the magnetic equator, the perpendicular and parallel intensities 
and are represented respectively by the solid and dashed lines. Note that the 
curves for £=1.0 MeV reproduce the values for a 0 and b 0 in Table 6 (cf. eq. 1 
and 3). The horizontal axis changes scale at L =5. 
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Figure 9. Integral omnidirectional flux for the Saturn proton model at three energy 
thresholds. The fluxes at the magnetic equator (B e /B=1.0) and at high magnetic 
latitude (B e /B =0) are represented respectively by the solid and dashed lines. The 
horizontal axis changes scale at L=5. 
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Figure 10. Integral intensity for the Saturn proton model at three sample values of L 
(numbers adjacent to the curves). On the left the dependences on pitch angle a 
are shown for energy thresholds of 1, 10, and 100 MeV (solid, dashed, and solid 
lines) at each L. On the right spectra are shown for (solid lines) and for /„ 
(dashed lines) at each L (/j_ and / . are equal for the isotropic distribution at 
L =2.67). 
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Figure 1 1 . Integral omnidirectional flux for the Saturn proton model at three sample values of 
L (numbers adjacent to the curves). On the left the dependences on local 
magnetic field strength B are shown for energy thresholds of 1 , 10, and 100 MeV 
(solid, dashed, and solid lines) at each L. On the right spectra are shown for the 
magnetic equator (B e /B = 1.0, solid lines) and for high magnetic latitude (B e /B= 0, 
dashed lines) at each L . 



92 


3 e /s 3e/3 


5217-90-029 


39 


A. R. Hoffman 


Figure 12. Coverage of the data in magnetic field strength B and threshold energy (see E a 
in Tables 3 and 4) for the electron (left panels) and proton (right panels) models 
at three values of L for each. The residuals are distinguished by values of A (cf. 
eq. 8) in the ranges <-0.6 (symbol ▼), -0.6 to -0.3 (symbol V), -0.3 to 0.3 (symbol 
o), 0.3 to 0.6 (symbol A), and >0.6 (symbol A). 
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APPENDIX VII. Review RFA Comments 


COMMENT #1 

Reviewer: Chris Paranicas 

Organization: Applied Physics Laboratory 
Topic: Energetic charged particle fluxes 

Concern: 

Updating the Saturn radiation model. 

Some Voyager channel fluxes have been corrected since the encounters. For example, PL04 is a 
total ion channel on Voyager but likely has been contaminated with electrons close to the planet. 
This means the previously reported PL04 rates are too high. 

Recommendation: 

Use the Cassini Low Energy Magnetospheric Measurement (LEMMS) measurements to improve 
the ion fluxes close to Saturn. 

As an additional action, it is worthwhile to compare the model predictions to the pure proton 
channel (channel 1) on the Voyagers near Saturn. 

COMMENT #2: 

Reviewer: W. Kent Tobiska 

Organization: Space Environment Technologies 
Topic: Standardization of Divine Model 

Concern: 

While many concerns and improvements to the Divine Model have been identified, it is likely 
that, due to programmatic constraints, this model will become widely used for engineering 
applications. This implies “standardization” of the model. Concerns are: 

1) How to ensure a reliable reference model for computational usage (published references, 
accessible code)? 

2) How are uncertainties (formal) and “weather” to be characterized? 

3) What is the model update strategy? 

Recommendation : 
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1) Use JPL code as a reference code once validated. Publish beyond JPL IOM in peer-reviewed 
literature. 

2) Clearly explain formal uncertainties in the model and how to account for space weather 
phenomena such as magnetospheric time variability and local time effects. 

3) Describe an update strategy — proposed or actual. 

COMMENT #3: 

Reviewer: Tom Armstrong 

Organization: Fundamental Technology, Inc. 

Topic: Reconciliation with Cassini SOI 

Concern: 

SATRAD should be compared with Cassini SOI observations of MIMI, CAPS, and 
Magnetometer. 

Recom mendati on : 

1) Obtain data (calibrated) and prepare graphs. 

2) Evaluate the results 

COMMENT #4: 

Reviewer: Tom Armstrong 

Organization: Fundamental Technology, Inc. 

Topic: Inclusion of PDS Archive Data 

Concern: 

Much more complete observations than used by Divine originally are available and should be 
used. Emphasize the 2-4 R s region. 

Recommendation: 

1) URL of VI and V2 Saturn PDS data available from Tom Armstrong. 

COMMENT #5: 

Reviewer: Tom Armstrong 

Organization. Fundamental Technology, Inc. 

Topic: Magnetic Field Models 
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Concern: 


1) Account for magnetopause currents. 

Recommendati on : 

1) Consider using Engle and Beard spherical harmonic expansions. 

COMMENT #6: 

Reviewer: Jack Connemey (Apr 22, 2004) 

Organization: NASA Goddard Space Flight Center 
Topic: Saturn Magnetic Field Models 

Concern: 

1) Does Saturn Z3 model in "Magnetic Fields of the Outer Planets" include external terms? 

Z3 is the best one can do for internal field. You're right, I didn't mention the external field 
coefficients, in the Z3 paper, they're not very useful in characterizing the external field for 
magnetospheric field line geometry. For that I'd direct you to a ring current model. 

2) The magnetic field goes to 0 around 6-10 R s and the field line tracing routine wanders off to 
infinity if it passes this region. 

The external field is something like 8-10 nT, as I recall, on average and on axis. You will 
get in trouble when the equatorial Z3 field in the theta direction approaches this value. But you 
ought to get well beyond 6 Rs before this happens (21,535/r**3), so if you have trouble inside of, 
say, 10 R s , something is wrong; look to external field coefficient normalization, or something 
that would alter the relative magnitude of intemal/extemal terms. 

3) Has any further work has been done with the Saturn field since Z3 and SPV; in particular are 
there any field models that use "ring currents" instead of the external harmonic terms? 

We did a ring current model early on, just like the Jovian magnetodisc model, and it works well 
to describe the field in the region where axisymmetry is OK - i.e., well within a magnetopause 
stand-off distance. You have code at JPL to do the calculation, just substitute the Saturn 
coefficients in the magnetodisc code. 

For further information, see the following references: 

"Zonal harmonic model of Saturn's magnetic field from Voyager 1 and 2 observations," J. E. P. 
Connemey, N. F. Ness, and M. H. Acuna, Nature, 298, 44 - 46, 1982. 

"Currents in Saturn's magnetosphere," J. E. P. Connemey, M. H. Acuna, and N. F. Ness, 
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J. Geophys. Res., 88, 8779 - 8789, 1983. 


If you want to code it up, see the Jupiter paper 

"Modeling the Jovian current sheet and inner magnetosphere," J. E. P. Connerney, M. H. Acuna, 
and N. F. Ness, J. Geophys. Res., 86, 8370 - 8384, 1981 . 


COMMENT #7: 

Reviewer: John F. Cooper ( 1 5 March 2005) 

Organization: NASA Goddard Space Flight Center 
Topic: Data concerns with the Divine Model 

Concern: 

My previous comments are now 1 5 years old and I have not spent much time revisiting these 
data since then, with a few exceptions as noted below. 

(1) 1 have long given up on the idea that LECP might be responding to a mono-energetic electron 
spectrum as suggested originally from the satellite filter model of Van Allen. My work on drift 
shell modeling for energetic electrons, to be reported at the Spring American Geophysical Union 
meeting, shows that such spectral features would be rapidly dispersed. The morphology of local 
time asymmetry is apparent in various Pioneer and Voyager electron data sets, and presence of 
these signatures in "nominal" proton or ion channels may suggest that these channels are 
responding to electron background within the inner magnetosphere. In general one should pay 
close attention to potentially anomalous morphology in counting rate data before such data are 
converted to fluxes .... 

(2) I remain concerned about potentially or actually anomalous responses of the Pioneer and 
Voyager particle instruments to the high-energy proton environment of the inner magnetosphere. 
This was very clear in the University of Chicago instrument L1NL2 > 0.5 MeV proton channel 
on Pioneer 1 1 and could be a concern elsewhere as possibly in the LECP high-energy electron 
channel > 10 MeV. Since MeV electrons are a critical issue for radiation hazards to spacecraft 
systems, you should certainly address this issue. 

(3) The LECP group gained much experience in correction of ion channels for electron 
background from Voyager measurements at Uranus and Neptune, and they have since said in 
publications that the earlier findings of heavy ions in the inner magnetosphere were erroneous 
due to this background. I note, however, that Uranus and Neptune did not have significant high 
energy proton components, so corrections inwards of Enceladus, where Voyager Cosmic Ray 
Subsystem last detected these protons, may still be uncertain. 

(4) It is crucial for the high-energy proton data to correctly map in B and L, since there is a 
northward dipole offset and this shows up very well in apparent radial offsets of inbound and 
outbound data if the mapping is not computed correctly. Assume that the high-energy proton 
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drift shells are longitudinally symmetric but that lower-energy electron and ion drift shells may 
not be. There can be both temporal and local time variations for the latter. 

I assume that you are setting the stage in this review for a model upgrade with all the new data 
from Cassini, although this mission did not have ALL the instruments really needed to fully 
characterize the high-energy particle environment of the inner magnetosphere. In particular, new 
data from a high-energy electron spectrometer would have been of great value to you now. So 
some uncertainties will remain. 


COMMENT #8: 

Reviewer: John F. Cooper (29 March 2005) 

Organization: NASA Goddard Space Flight Center 
Topic: Data concerns with the Divine Model 

Concern: 

My main comment re Chicago Pioneer 1 1 data would be that the "0.5-10 MeV proton" 
(L1NL2) and "7-17 MeV electron" (D5ND67) channels from the Simpson instrument appear by 
my previous analyses to have mixed nominal and high-energy (>30 MeV) proton responses in 
the inner magnetosphere inwards of Mimas and outwards from the outer edge of the A ring, so 
these channels should not be given much weight in that region. I would also be suspicious of the 
> 10 MeV electron rate from Voyager LECP in that region for the same reason. In this region, I 
would attribute inbound-outbound differences from the more equatorial Pioneer 1 1 data to local 
time asymmetry associated with drift resonance effects rather than to temporal variations. If the 
local time asymmetry is driven by a dusk-to-dawn electric field, then one would expect higher 
energetic electron intensities in the dusk-centered hemisphere, as compared to strong 
longitudinal symmetry for the high-energy protons as mapped in the offset-dipole system. 

Finally, the region under the rings is probably well represented by the time-averaged 
proton, electron, neutron, and gamma ray spectra in the Cooper et al. (1985) paper from the 1979 
Pioneer 1 1 data. 

If this is possible, I certainly would like to request access to the Pioneer-Voyager Saturn 
model in support of my joint work with Ed Sittler (GSFC) and Bob Johnson (Univeristy of 
Virginia) of the Cassini Plasma Spectrometer (CAPS) team on magnetospheric interactions with 
rings and moons as sources of plasma and neutral molecules for the magnetosphere. To the 
extent that this work also involves work with Cassini plasma and energetic particle data, I would 
be happy to contribute to a new Cassini-based model. 
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APPENDIX VIII. Acronyms and Abbre\nations 


APL 

CTS 

ESA 

FT 

GIRE 

GSFC 

JHU 

KU 

LECP 

MIMI 

MeV 

NAIF 

PDS 

RFA 

R j 

R s 

SATRAD 

SET 

SOI 

UT 


Applied Physics Laboratory [the Johns Hopkins University] 
counts per second 
European Space Agency 
Fundamental Technologies , Inc. 

Galileo Interim Radiation Electron 

NASA Goddard Space Flight Center 

Johns Hopkins University 

University of Kansas 

Low Energy Charged Particle Instrument 

Magnetospheric Imaging Instrument 

million electron volts [unit of energy] 

(NASA) Navigation and Ancillary Information Facility 
Planetary Data System 
Request for Action [form] 

Jupiter radius [71,400 km] 

Saturn radius [60,400 km] 

Saturn Radiation Model 
Space Environment Technologies 
Saturn Orbit Insertion 
Universal Time 
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